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INTRODUCTION 

The SOFIA Aircraft System (AS) Phase I Study was a broad-based 
study which addressed itself to satisfying technical and 
programmatic requirements by drawing from existing technology and 
applying cost-efficient commercial approaches to the aircraft 
modification. 

In this SOFIA AS Phase II Study, five (5) critical areas of the 
aircraft were selected for more detailed investigation. These are; 

Forward Pressure Bulkhead 
Aft Bulkhead 

"Free” Shell to Bulkhead Interface 
Shell Cut-out Surround Structure 
Flooring 

The in-depth investigation of these areas upon a Finite Element 
Model (FEM) , with a Fine Grid Model in areas of particular 
interest. The FEM code used is called "STRAP” and was developed by 
the engineering firm Rasmussen and Associates. STRAP is NASTRAN 
compatible to within 1%. The loads applied to the model were 
approximated from known 747 envelope conditions. The areas of 
investigation, and a section through the fuselage is shown on pages 
ii and iii. 

The thrust of this investigation was to develop the design concepts 
conceived under phase I to the point where detail design could be 
undertaken with a high confidence level. 
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(SURROUND STRUCTURE) 


PRINCIPAL AREAS OF STUDY AFT CAVITY LOCATION 
SOFIA PHASE II 
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VIEW SHOWING SOFIA TELESCOPE LOCATION 
IN 747-100 SIZE AIRCRAFT 


PC005-045 



FORWARD PRESSURE BULKHEAD 
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This bulkhead forms the aft boundary of the pressurized portion of 
the fuselage and provides primary support for the Telescope 
assembly. It also forms the forward wall of the telescope cavity. 

The bulkhead is constructed of twelve built-up 20 inch deep beams, 
four oriented vertically, four horizontally, and four on the 
diagonal. At the intersections the beams are spliced using cap 
slice plates and "bathtub” fittings. These beams support a flat 
extruded pressure web with bulb risers facing aft into the cavity. 
An extruded "Tee" section cap provides the web to fuselage skin 
attach member. 

The bulkhead is stiffness critical as a result of a criteria of 0.5 
inch deflection at the telescope mount points. The member sizes 
required to achieve this stiffness result in low stress levels for 
all the applied load conditions. AFT BULKHEAD 

This bulkhead is of similar construction to the forward bulkhead 
although of less massive proportions as it does not have to resist 
the cabin pressure loads and carry the weight of the telescope. It 
does however carrry the telescope brake, a spring/hydraulically 
operated device, and so is "tuned" to the stiffness of the forward 
bulkhead to minimize brake stroke differences between the 
pressurized and unpressurized conditions. The bulkhead also 
contains a non-structural door for ground access to the telescope 
cavity from the aft fuselage. 


BULKHEAD TO SHELL INTERFACE 

The attachment of the forward bulkhead attach to the fuselage shell 
is a critical fatigue joint. The flat bulkhead rotates aft at the 
shell intersection and constrains the shell against expansion while 
forward of the bulkhead the shell expands normally under cabin 
pressure. 

Intercostals of tailored stiffness to the shell at each twenty inch 
deep beam to provide a gentle transition from the bulkhead 
constrained shell to the "free" shell. Due to the critical nature 
of this joint a typical intercostal was modeled down to the 
fastener level. The detailed loads extracted from this model run 
were then used to calculate Stress Severity Factors (SSF) which, in 
turn, were used to determine a fatigue life for this typical joint. 
The intercostal was sized to give a predicted factored life of 
greater than 100000 cabin pressurization cycles. 
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CUTOUT SURROUND STRUCTURE 

The purpose of the cutout surround structure is to restore the 
strength and stiffness lost due to the shell cutout. This is done 
by framing the opening with longerons and using the bulkhead caps 
and fuselage frames as the other framing members. Skin doublers 
and triplers are added to accommodate the increased shear flows and 
provide added stiffness and stability at the cutout edge. Due to 
the anticipated configuration of the NASA designed cavity door the 
primary longerons (the ones immediately adjacent to the cutout) are 
external to the fuselage and aerodynamically faired, the secondary 
longerons are internal. As the longerons are essentially 
unsupported between bulkheads, they require high moments of 
inertia, particularly the primary ones. 


FLOORING 

with the exception of a maintenance walkway the aircraft main deck 
floor is removed between the bulkheads to permit installation of 
the telescope. A walking surface is also added to the cargo floor. 

Forward of the pressure bulkhead and aft of FS 1600 the floor is 
cut away between BL 60 L & R to permit installation of the 
telescope and provide comfortable access to the telescope 
instrument package. The retained portion of the floor has shear 
webs and fore/ aft beams added to enable it to act as horizontal 
beams to hold the shell in the same shape as when it is tied across 
by floor beams. 

The floor is also cut back from the bulkheads to permit their 
deflection without imposing loads on the floor structure. 


CONCLUSIONS 
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Structurally the SOFIA program is very "doable”. It can be 
achieved using standard aircraft materials and processes and does 
hot require exotic materials to achieve its goal. As can be seen 
from the preceding section, stiffness is as much a driving force as 
strength. This does tend to make the structure rather heavy, but 
as the telescope assembly, support structure, and instrumentation 
package comprise most of the effective payload weight, this is not 
a significant problem for this aircraft. 

Structurally the most problematic aspects of this program are: 

o Boeing unwillingness to furnish data, and hence the need to 
use a "back-door" approach, (a frequently used option in the 
after-market modification field) . 

Comment: This is not difficult to overcome although it will 

increase engineering costs. 

o Fatigue damage at the attachment of the flexible shell to the 
rigid bulkhead. 

Comment: This can be addressed through attention to detail 

design and adequate component testing. 
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PROPOSED PHASE II STUDY TOPICS 

As stated in the Conclusions section, the most problematic 
structural item is the bulkhead to shell joint. In light of this, 
a Phase III Study could address this issue with a detail design 
study coupled with a full scale fatigue test of the joint. 
Stringer fitting feed-throughs should be incorporated in the test 
item to provide a comprehensive fatigue assessment of the system. 

A second and related study would be a detail design and fatigue 
test of the pressure bulkhead beam intersections. 

LAS recommends these subject areas be included in a Phase III SOFIA 
Aircraft System study. 
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ABSTRACT 

This report furnishes the results of the Phase II study contract 
for the Stratospheric Observatory for Infrared Astronomy (SOFIA) as 
installed in a 747-100 aircraft. 

The SOFIA mission requires the airborne telescope to be Open to the 
atmosphere during observation periods and isolated from ambient 
conditions during pre-flight periods for cooling. These 
constraints require a large cut-out in the fuselage for viewing and 
bulkheads for support and isolation, both forward and aft, of the 
telescope. A closeout door is also required to isolate the cavity 
during the pre-cooling period. The design of this item has been 
retained by NASA. 

This study is limited to five specific areas as defined by the 
statement of work, these being: 

i ) Forward Bulkhead 

ii) Aft Bulkhead 

iii) Shell Cutout Surround Structure 

iv) Forward Bulkhead to Shell Transition 
V) Floor Structure 

The criteria, assumptions, methodology, and results are described 
in the following sections and appendices. 
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1.0 INTRODUCTION 


The SOFIA installation investigated in this report uses a 747-100 
aircraft as its platform; this is at variance with the Phase I 
study which used a 747 SP as the transport system (reference LAS 
report R3604001) . This change was made to minimize cost on this 
limited-budget study. Lockheed had access to a previously 
developed model and loads which had been derived from, and FAA 
approved for, a 747-100 modification program. Development of loads 
for a 747SP SOFIA aircraft would have exceeded the available 
budget, even for a very simplified set of loading conditions. 
Also, should NASA elect to use a 747SP as the system platform, the 
bulk of the results contained in this report would still apply in 
principle; they need only be modified by the results of the detail 
design process as they would (although to a lesser degree) for a 
747-100 platform. The basic 747-100 aft fuselage geometry is shown 
on page 1.2. 

The fundamental SOFIA configuration is based on input from NASA and 
Lockheed Missiles and Space Company SOFIA Final Review 
LMSC/P016023 . The required fuselage shell opening was defined by 
NASA and the Telescope mass properties were obtained from the LMSC 
report . 

The general arrangement of the installation is shown on page 1.3. 
The bulkhead to bulkhead distance is 180 inches, the shell cut-out 
fore/aft distance is 160 inches and it is 89 degrees wide 
circumferentially. With these parameters the aft face of the 
forward bulkhead was set at FS 1780 and the forward face of the aft 
bulkhead at FS 1960. The center rotation of the telescope is at FS 
1770, WL 231 (the origin of the upper fuselage's 117.75 inch 
radius) and BL 0. Between the bulkheads, the main deck floor is 
removed with the exception of a maintenance walkway, to permit 
installation of the telescope. Forward of the forward bulkhead and 
aft of FS 1600, the floor is cut away between BL 60 L & R to permit 
installation of the telescope mounting and operating mechanisms and 
to permit access to the instrument package. 

The aircraft retains its normal pressurized cabin forward of FS 
1780. The area between the bulkheads is protected from the 
external environment by a moveable door and orifice (NASA design) , 
and is open to atmosphere during observation periods. The fuselage 
aft of the aft bulkhead is vented to atmosphere. 


FRONT SPAR BHD 
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INTRODUCTION (cont'd) 
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2.0 CRITERIA AND ASSUMPTIONS 

The criteria and assumptions on which the SOFIA Phase II study was 
based are delineated in this section. 

2.1 ALLOWABLE STRESS LEVELS 

Several magnitudes of allowable stress levels apply to the SOFIA 
installation depending upon the function of the part, and upon the 
applied loading. 

2.1.1 FATIGUE CRITICAL PARTS 

The allowable stress level for fatigue critical parts are gross 
area tension stresses assuming a joint efficiency of 75% i.e. no 
more than 25% of the tension area removed by fastener installation. 
The stress levels are based upon Lockheed developed constant life 
curves for standardized aluminum. 

Inspectable installations which are subject to flight loads with or 
without pressure loads: allowable gross area ultimate stress level 
is 45 ksi. As no fatigue spectrum is available to scientifically 
derive a value for this condition it is based upon past experience 
with wide-body commercial aircraft. 


Uninspectable areas of installations which are subject to flight 
loads with or without pressure loads: allowable gross area 
ultimate stress is 30 ksi. As with the 45 ksi allowable above, 
this is based upon past experience. 

Built-up structure subject to cabin pressurization only loads: 
allowable gross area tension stress is 14 ksi, limit. This 
allowable is based on a fatigue quality index of 4.5 and factored 
life of approximately 50000 cabin pressurization cycles. The 
allowable stress is low as it is assumed that full cabin 
pressurization (9.4 psi) is achieved every flight. The allowable 
was determined using the Constant Life Curves for standardized 
aluminum from the Lockheed Structural Life-Assurance Manual 
(SLM)and the results are. shown on page 2.2. 

The outstanding flanges of the bulkhead pressure panels which are 
subject to cabin pressurization loads: allowable gross area 
tension stress 30 ksi, limit. This allowable is based on a fatigue 
quality index of 1.5 and hence, no attachments are permitted in the 
riser beyond the riser mid-height. This was determined from the 
SLM in a similar manner to the built-up structure, and the results 
are shown on page 2 . 3 
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2.0 CRITERIA AND ASSUMPTIONS (cont'd) 

2.1 ALLOWABLE STRESS LEVELS (cont'd) 

2.1.2 Non-Fatigue Critical Parts 

For non-fatigue critical structure the allowable stresses are those 
given in MIL-HDBK-5F or other FAA approved sources. Yield 
allowables shall not be exceeded under limit load conditions and 
ultimate allowables may not be exceeded under ultimate conditions 
(i.e. 1.5 X limit load or Emergency Landing Conditions per FAR 
25.561). Plastic bending may be used for ultimate conditions 
providing this does not conflict with the yield criteria at limit 
load. 

For crash conditions the structure may suffer permanent deformation 
but shall not fail. 

2 . 2 STIFFNESS REQUIREMENTS 

The overall vertical and lateral bending stiffness of the SOFIA 
modified aircraft shall not be less than the stiffness of the 
unmodified aircraft. The torsional stiffness shall be as high as 
possible, within reason, with equivalence to the unmodified 
aircraft being the goal. 

The forward bulkhead which supports the telescope shall deflect aft 
no more than 0.5 inches under limit loads at the telescope mounting 
points. The aft bulkhead stiffness shall be tuned to that of the 
forward bulkhead such as to minimize deflection effects from the 
pressurized to the unpressurized conditions on the telescope brake 
system. The brake is assumed to be actuated by a spring with 
hydraulic brake releases. The brake is assumed to be engaged under 
hydraulic failure. 

2.3 MISCELLANEOUS LOAD DERIVATION AND ASSUMPTIONS 

The telescope mass properties are taken from the LMSC Final Review 
for the telescope design (Report No. LMSC/P016023 dated 6/18/93) 
and are repeated here for convenience. 

Telescope System Weight = 12861 kg 

Telescope Inertia (kg-m2) ^ 

6051 77 

I = 77 39733 

-2595 58 


-2595 

58 

39485 


From ref 3 Preliminary FEM Report 

Max Yaw Acceleration = 0.212 rad/sec2 
Max Pitch Acceleration = 0.363 rad/sec2 
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2.0 CRITERIA AND ASSUMPTIONS (cont'd) 

2.3 MISCELLANEOUS LOADS DERIVATION AND ASSUMPTIONS (cont'd) 
Brake Load required to restrain telescope 


As the above loads are preliminary, conservatively use a rotational 
acceleration of 0.5 rad/sec2 for brake sizing. 


Max inertia = 39733 kg-m2 
Max yawing moment 
Bearing to Brake Face arm 
Load normal to brake axis 
Friction coefficient 
Axial brake force 


= 29205 lb-sec2-ft 
= 14602.8 Ibft 
= 155 in 
= 1130 lb 
= 0.23 (assumed) 

= 4913 lb (limit) 

= 7370 lb (ult) - use 7500 


For design study purposes the ultimate brake load is 

7500 lb axial 
1725 lb radial 


An arbitrary 1.0 psi limit load is assumed to be acting on the aft 
bulkhead (1.5 psi ult) 
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3.0 FINITE ELEMENT MODEL 

3.1 BASIC MODEL DESCRIPTION 

SOFIA is modeled using the "STRAP” finite element modelling code. 
STRAP is a NASTRAN compatible code developed by Rasmussen and 
Associates, Burbank, California. The modelling was achieved by 
first modelling a standard 747-100 aircraft as described in the 
following sections. 

3.1.1 GEOMETRY 

The fuselage of the Boeing 747-100 aircraft aft of station 1480 is 
a nearly cylindrical monocogue with a modest upward sweep and 
contraction aft at approximately Station 1780. The Station 1480 
Bulkhead is located just forward of the model, and the wheel well 
is located just ahead of the bulkhead. The wheel-well side lobes 
(left and right) and keelson attach to the bulkhead. The presence 
of the discontinuity in the monocogue structure, represented by the 
wheel well and the loads introduced by the side lobes and keelson 
cause a redistribution load that is noticeable up to about 300 
inches aft of the bulkhead. For this reason, an adeguate 
representation of the structure ahead of Station 1480 is necessary. 

3.1.2 MONOCOQUE 

The skin is supported by 51 stringers per side and the stringers 
are supported by frames with a 20 inch spacing throughout. 
Stringers 41 and 43 terminate at Frame 2120. Stringers 27 and 45 
terminate at Frame 2140. Stringer 39 terminates at Frame 2160. 
The finite element model contains grid points at all stringers to 
frame connecting points. All frames and all stringers are modeled 
without lumping, except stringer 32A which extends (within the 
model) from Station 1480 to Station 1520 and is slightly offset 
from Stringer 33. This stringer is lumped with Stringer 33. The 
skin geometry is shown on Page 3.4. The stringer locations and 
the skin and stringer properties have been taken from Boeing data. 
These properties closely match those used by Boeing in the 
monocogue stress analysis, however, the box-beam analysis performed 
by Boeing will reduce the effectiveness of skin segments in 
compression, and will introduce compressive loads in stringers due 
to panels in diagonal tension. To capture these phenomena in a 
finite element analysis, the stiffness properties would have to be 
modified depending on the load condition; i.e. separate upbending, 
downbending, lateral bending, and pressure models would have to be 
used. These refinements are not reguired for a preliminary study 
and have not been implemented. Internal loads at Station 1740 for 
the FE model without entry doors compare reasonably well with 
Boeing values (see Page 3.19 and 3.20). 
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3.0 FINITE ELEMENT MODEL (cont'd) 

3.1.3 TYPICAL FRAME 

The frame dimensions follow the skin mold line as defined for the 
monocogue geometry. A cross section at Station 1740 is given on 
Page 3.5. Frames are located 20 inches apart from Station 1480 
through Station 2180 at both end of the aft fuselage model. The 
stiffness properties of the frames are those of the representative 
frame at Station 1600 from Boeing data. 

3.1.4 FLOORS 

The aft fuselage has two floors, the main deck at WL 199.3, and the 
cargo deck at WL 118.6. The main deck is shown on Page 3.6 and its 
aft end on Page 3.7. The main deck is supported by transverse 
floor beams at each frame. The transverse floor beams are 
supported by posts (left and right) at BL 96.0. Seat tracks are 
located at BL 11.33, 33.99, 75.92, and 96.0, with auxiliary seat 
tracks at BL 67.01 and 89.67 aft of Station 2080. The seat tracks 
at BL 33.99 and are mounted on top of a longitudinal floor beam. 
This floor beam is continuous forward of Bulkhead 1480, whereas the 
other seat tracks are structurally interrupted at Bulkhead 1480. 
The cargo deck is supported directly by the inside cord of the 
frame, which is a deep beam between BL 60L and 60R. The cargo deck 
slopes up at the forward end of the bulk cargo compartment (aft of 
Station 1920). (Reference Page 3.8) 

3.1.5 ENTRY DOORS 

The aft fuselage model contains the entry doors (left and right) 
between Frame 1660 and 1720 and the Cl Lower Hold Cargo Door 
between Frame 1800 and Frame 1920 on the right hand side. The 
lower Hold Bulk Cargo Door aft of Station 1980 on the right hand 
side has been omitted. Entry doors have been omitted and openings 
have been modeled as "skinned over" with the original skin, 
stringer, and frame properties. The surround structure of each 
entry door consists of upper main, upper auxiliary, lower main and 
lower auxiliary skills, frame reinforcements, and skin gage 
increases. Properties have been obtained from Boeing data. 
(Reference Page 3.9 & 3.10) 

3.1.6 SUBSTRUCTURE FORWARD OF STATION 1480 

The fuselage and wing forward of Station 1480 has been modeled as 
a substructure (or super-element) . This substructure represents 
the Station 1480 Bulkhead, the stub wings and the fuselage center 
section in the analysis of the aft fuselage. The general lay-out 
of the FE model ahead of Station 1480 is shown in this report. 
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3.0 FINITE ELEMENT MODEL (cont'd) 

3.1.6 SUBSTRUCTURE FORWARD OF STATION 1480 (cont'd) 


The substructure has been derived from a Boeing 747 wing to 
fuselage intersection model. The accuracy of this model is not 
guaranteed for the individual stiffness elements. However, the 
overall stiffness properties of this model are believed to be 
sufficiently accurate for the purpose of constraining the current 
aft fuselage model at Bulkhead 1480 by means of the constraint 
points provided in the forward fuselage and at the stub-wing tips. 
A view of the model is shown on Page 3.11. 

The baseline model from which SOFIA was developed is shown 
pictorially on the following pages as listed on the chart below: 


PAGE 

DESCRIPTION 

3.4 

Side view of Fuselage Shell FS1480 to 2180 

3.5 

Cross section of a Typical fuselage frame (FS1740) 

3.6 

Top view of main deck floor (FS2040 to FS2180) 

3.7 

General view of main deck floor (right hand side) 
FS1480 to FS2180 

3.8 

General view of cargo floor (right hand side) 
FS1480 to FS2180 

3.9 

Skinned over cargo door and framing reinforcement 

3.10 

Skinned over passenger entry door and framing 
reinforcement 


3.11 


General view of sub-structure forward of FS1480 
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3.0 FINITE ELEMENT MODEL (cont'd) 

3.2.0 LOADS 

Preliminary design loads for the aft fuselage model have been 
developed from envelope loads in Boeing data. The shear moment and 
torque distribution of the load cases applied to the model are 
shown on Pages 3 . 17 and 3.18. 

3.2.1 LOAD CASE SELECTION 

Three basic load cases have been developed; a vertical down bending 
envelope case, a pure lateral bending envelope case, and a pressure 
case for pressure relief valve pressure difference of 9.4 psi. The 
bending envelope cases have shears and moments introduced at each 
frame station of the model to match the envelopes from Boeing data 
in both shear and moment just aft of each frame. The three basic 
load cases are combined into five additional load cases. First, 
the vertical load case is scaled to obtain an approximate 1.5-g 
vertical load distribution, representing level flight. The scaling 
factor is 1.5/3.75, because the vertical load envelope is 
determined by maneuver conditions having a 3.75 load factor over 
most of the aft fuselage length. The scaled vertical bending 
envelope is combined with the lateral load case to obtain the 
fourth load case; lateral load combined with level-flight. A fifth 
load case is obtained by scaling the pressure conditions to 2.0 
factored pressure. The three remaining conditions are obtained by 
combining 1.5 factored pressure with the two lateral and the 
vertical bending cases. 


Condition 

Description 

Vertical 

Downbending Envelope 

Lateral 

Lateral Bending Envelope 

IP 

9.4 Psi Cabin Pressure (P) 

Ig FLT 

1.5g Level Flight 

LateNzo 

1.5g Level Flight & Lateral Bending 

2P 

18.8 Psi Cabin Pressure 

Vertical + P 

1.5P + Downbending Envelope 

Lateral + P 

1.5P + 1.5g Level Flight & Lateral 
Bending 

LateNzo + P 

1.5P + Lateral Bending 
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3.0 FINITE ELEMENT MODEL (cont'd) 

3.2.1 LOAD CASE SELECTION (cont'd) 

Internal loads at Station 1740 are compared to Boeing Box Beam 
analysis results on Page 3.18 and 3.19 for the vertical envelope 
condition "Vertical" and the lateral envelope conditions "Lateral" 
and LateNzo", representing lateral loads with and without vertical 
loads, respectively. The internal loads for the Vertical Condition 
show a fairly good match, both in shear flow and in longitudinal 
stress, except that longitudinal stress exceeds Boeing values in 
the top (in tension) and is somewhat less in the bottom (in 
compression) of the monocoque. These difference are considered a 
result of differences between Box Beam and FE analyses, as 
mentioned in Section 3.1.1. The Boeing lateral condition is 
bracketed in shear between the lateral conditions with and without 
1-g vertical loads. The longitudinal stress for the lateral 
condition with vertical shear is slightly exceeded by the stress 
for the Boeing condition (in compression) around Stringer No. 30. 
The Boeing box beam analysis is simply a unit beam solution with 
adjustments made to accommodate cutout effects. 

3.2.2 LOAD CASE APPLICATION TO MODEL 

The incremental shears, bending moments, and torsion are applied to 
reference points at each frame in the plane of symmetry at W.L. 
225. These points are connected to the fuselage monocoque 
structure using a splining element (RSPLINE) that introduces the 
loads without modifying the stiffness of the model. This element 
employs a factored least-squares technique to distribute loads into 
each frame. Shears are distributed uniformally in pounds/ inch over 
the circumference, and moments are distributed as a linear 
pounds/ inch distribution about a neutral axis. Torsion is 
introduced as a nearly uniform T/2A shear flow. Most of the 
incremental loads introduced at the frames consist of inertia 
loads, as airloads on the fuselage are relatively small. 
Therefore, the true application pattern of the inertia loads is 
less smooth than the one employed here. In both cases, smooth or 
discrete, the stiffness properties of the frame and adjacent skin 
cause the applied incremental loads to be redistributed in 
accordance with the compliance of the structure within a frame 
spacing length. In the Boeing box beam analysis, incremental loads 
are also applied in a smooth distribution. 

For variations of the current model, the splining can easily be 
modified by removing or adding load points in the definition. The 
distribution of load can be influenced by changing factors. For 
example, large masses, such as the SOFIA telescope, can be 
represented by interior points with a relatively high factor (to 
attract load) . 
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3.0 FINITE ELEMENT MODEL (cont'd) 

3.2.2 LOAD CASE APPLICATION TO MODEL (cont'd) 

The benefit of the splining method is that load can be distributed 
into complex models with little effort. However, an acceptable 
load distribution for (preliminary) design can only be obtained 
with a sufficient number of points, distributed vertically and 
horizontally over the frame. A point set without vertical spread, 
for example, cannot be used to introduce vertical bending. The 
load distribution that is created can be examined using grid point 
force analysis. 

Pressure loads are applied to the model using a distributed 
pressure over each panel. At the model ends, loads are distributed 
to triangular areas: two support points along the edge of the 
monocoque, and the third point a centrally located reference point. 
The Siam of the load applied to the reference point is transferred 
back to the edge using the splining method described above. Again, 
the total pressure load is area-integrated in the finite element 
calculation. Subsequently, this load is distributed in a uniform 
manner over the edge of the monocoque by the spline. At the 
forward end, pressure loads are applied at the wheel-well pressure 
bulkhead only. The balance of the pressure load, equal to the 
pressure times the open area above the wheel well at Station 1480, 
is reacted by the substructure forward of that station. 

3.2.3 AFT FUSELAGE ACCELERATION 

Critical SOFIA design vertical and lateral accelerations are 
obtained for light-weight gust conditions and a two-wheel level 
landing condition. The actual SOFIA A/C weight and c.g. location 
have not yet been determined. Therefore, an Operating Empty Weight 
(OEW) of 327,000 lbs. at 34.9% from Boeing data has been 
incremented by 33,000 lbs. of weight in the forward fuselage to 
bring the c.g. to the forward limit of 14% MAC. Yaw and pitch 
inertia have been estimated based on the length of the fuselage, 
the OEW and the fuel and payload weights are listed on Page 3.15. 
Actual SOFIA weights will be higher than those in the table and 
will have a more aft c.g. due to the telescope installation weight 
that is added to the aft fuselage. 

Therefore, actual accelerations will be less than those calculated 
in this report. Accelerations are documented on page 3.15. The 
gust accelerations have been calculated using FAR 25.341(b) 3 and 
25.1351(b) for gust at Vc and the aerodynamic derivatives from the 
table on Page 3.16. 
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3.0 FINITE ELEMENT MODEL (cont'd) 

3.2.3 AFT FUSELAGE ACCELERATION (cont'd) 


747-100 Weight Conditions 


Cond. 

GW 

kips 

c.g. 

MAC 

lYY 

in=lb-sec2 

Izz 

in-lb-sec2 

000 

327 

34.9% 

.258E+09 

.374E+09 

001 

360 

13.0% 

.305E+09 

.422E+09 

002 

585 

13.0 

.480E+09 

.611E+09 


747-100 LIMIT Acceleration Design Conditions 
(NY, 'Pddot apply to Lateral Gust, Nz, i)ddot to others. 
Gust Conditions at Vc; 375 KEAS at 20kFT, or M 0.92 at 
40 kFT) 


Cond. 

Type 

GW 

C.g. 

Alt. 

Ny/Nz 




kips 

%MAC 

KFT 

LIMIT 

rad/ sec2 

GLOl . 1 

Lat. Gust 

360 

0.13 

20 

0.51 

0.212 

GL02 . 1 

Lat. Gust 

360 

0.13 

40 

0.24 

0.100 

GNOl.l 

-Gust 

360 

0.13 

20 

-1.02 

0.363 

GN02.1 

-Gust 

360 

0.13 

40 

-0.16 

0.201 

GPOl.l 

+Gust 

360 

0.13 

20 

3.02 

-0.363 

GP02 . 1 

+Gust 

360 

0.13 

40 

2.16 

-0.201 

LL01.2 

2 WLL 

585 

0.13 

S.L. 

2.15 

-0.190 


2 WLL Two-wheel Level Landing at Maximum Vertical 
Reaction 
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3.0 FINITE ELEMENT MODEL (cont'd) 

3.2.3 AFT FUSELAGE ACCELERATION (contM) 

747-100 Geometric Data 


Mean Aerodynamic Chord 

MAC 

327.8 

; in 

Leading Edge of MAC 

LEMAC 

1258 

in 

Wing Reference Area 

S 

5500 

ft-sq 

Wing Span 

B 

2348 

in 

Length of Fuselage 

LFUS 

2702 

in 

Fuselage Station of 
Horizontal Stabilizer 
25% MAC 

Xtail 

2495 

in 

Distance from 25% Wing 
MAC to 25% Horizontal 
Stabilizer MAC 

Ltail 

1254. 

05 in 

Horizontal Stabilizer 
Reference Area 

Sht 

1470. 

7 ft-sq 

Wing Aspect Ratio 

A 

6.96 

Horizontal Stabilizer 
Aspect Ratio 

Aht 

3.6 
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3.0 FINITE ELEMENT MODEL (cont'd) 

3.2.3 AFT FUSELAGE ACCELERATION (cont'd) 


747-100 Aerodynamic Derivatives 


Altitude 

Alt. 

(kft) 

20 

40 

Velocity 

Ve 

(keas) 

375 

262 

Mach Number 

M 

0.836 

0.92 

Lift Derivative 
Tail Off 

CL alpha 
(1/deg) 

0.0736 

0.0855 

Pitching Moment 
Derivative Tail Off 

Cm alpha 
(1/deg) 

0.0150 

0.0084 

Rate of Change of Down Wash 
with (respect to alpha) 

deda 

0.3857 

0.4481 

Lift Derivative Tail On 

CL alpha 
(1/deg) 

0.0852 

0.0959 

Pitching Moment Derivative - 
Tail On 

CM alpha 
(1/deg) 

0.0293 

-0.0314 

Side Force Derivative 

CYbeta 

(1/rad) 

-1.0258 

-1.0258 

Directional Stability 
(Weather Cocking Derivative) 

Cnbeta 

(1/rad) 

0.1964 

0.1964 


Where ALPHA is angle of attack 
BETA is angle of yaw 


Acceleration Due to Gravity 


g 


32.2 ft/sec-sq. 
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3.0 FINITE ELEMENT MODEL (cont'd) 

3 . 3 ENVELOPE LOADS 

3.3.1 Aft Fuselage Vertical Envelope Loads 



PC005-058 
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3.0 FINITE ELEMENT MODEL (cont'd) 

3.3 ENVELOPE LOADS (cont'd) 

3.3.1 Aft Fuselage Lateral Envelope Loads 
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3.0 FINITE ELEMENT MODEL (cont'd) 

3 . 4 OUTPUT COMPARISON 

3.4.1 Shear Flow Comparison - Aft Fuselage without Doors FS 1740 



STRINGER # 


P0005-061 
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3.0 FINITE ELEMENT MODEL (cont'd) 

3.4 OUTPUT COMPARISON (cont'd) 

3.4.1 Longitudinal Stress Comparison - Aft fuselage without Doors 
FS 1740 



STRINGER # 


PC006-060 
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3.0 FINITE ELEMENT MODEL (cont'd) 

3.5 SOFIA MODEL 

The basic 747-100 Finite Element Model described in the preceding 
sections was modified to reflect the SOFIA configuration. 

The stringers were eliminated at the shell cutout and the skin 
panel thicknesses increased to reflect the surround doublers. The 
stringer elements in the area of the primary and secondary 
longerons were increased in area and moment of inertia and also had 
offsets incorporated to simulate the longeron centroid shifts. In 
the initial modeling the longerons were discretely modeled with 
individual elements for caps, webs, etc., but with evolving designs 
this was found to be too time consuming considering the budgetary 
constraints, and so the lumped approach was taken. 

The forward bulkhead was discretely modeled throughout the study as 
it was felt to be important that the loads at the bulkhead to shell 
transition be understood. The transition intercostal was modeled 
down to the fastener level at one location while at all other 
locations they were incorporated as beam elements only. The 
discretely modeled intercostal was run for the 9.4 psi condition 
only, as this is the condition on which the life estimate is based. 

The aft bulkhead was modeled using beam elements for the bulkhead 
beams as opposed to discretely modeled elements. 

The aft fuselage model is shown on page 3.23. 





Page 4 . 1 
Model 747 
Report R54045007 

4.0 STRUCTURAL DESCRIPTION AND ANALYSIS RESULTS 

4 . 1 FORWARD BULKHEAD 

The forward cavity bulkhead is centered at FS 1770 and performs 
three major functions: 

i) It forms the aft boundary of the pressurized portion of the 
fuselage. 

ii) It supports the mass of the telescope and associated 
equipment. 

iii) It acts as the primary forward framing member for the 
observation cut-out. 

The bulkhead and support beams are shown on sheets 2 , 3 , 4 , & 8 of 
Drawing B54 114000. 

The bulkhead is a flat pressure web supported by 20 inch deep beams 
with the web located at FS 1780 and the support beams extending 
forward to FS 1760. The bulkhead contains a 100 x 100 inch cut-out 
to accommodate the telescope support bearing and pressure seal. 

The pressure web is fabricated from 20 inch wide extruded panels 
with a maximum length of 20 inches between supports. These panels 
have 0.94 inch high bulb risers at approx 1.5 inch 0/C (Lockheed 
extrusion no. LS5150) . The risers face aft into the cavity and are 
locally machined away to accommodate the cap splice plates and 
pressure plenum. The basic LS5150 extruded panel is 7075-T6 but as 
this design will have extensive machining to accommodate the splice 
plates and plenum seal, the extrusion should be in the T73 over- 
aged condition to minimize exfoliation and stress corrosion 
problems . 

The main beams are 20 inches deep and are of conventional cap, web, 
and stiffener construction. As the bulkhead is stiffness critical 
the shear web is designed to be non-buckling under normal operating 
conditions. This results in the shear webs outside the circle 
inscribed by the pressure plenum seal being 0.08 inch thick and 
those inside the circle being 0.063 thick. In addition, where the 
bulkhead to shell transition intercostals occur, the webs are 
locally doubled up to present a flat surface to the transition cap 
members. The caps are extruded '‘Tee” sections tapering in 
thickness toward the shell. The caps are spliced at the joints 
using external plates and internal bathtub fittings and bolts. 
This is true for both the forward and aft caps. 
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4.0 STRUCTURAL DESCRIPTION AND ANALYSIS RESULTS (cont'd) 

4.1 FORWARD BULKHEAD (cont'd) 

The recommended materials for this bulkhead are: 


Pressure Panel 
Shear Beam Webs 
Shear Beam Caps 
Splice Plates 
Splice Fittings 


LS5150 aged to T73 

7075 T6 Clad 

7075 T73 Extrusion 

7075 T73 Plate 

7075 T73 Hand Forged Billet 


The results of the FEA model runs are summarized on the following 
pages as shown in the chart below: 


Page 

Description 

4.3 

Beam Web Shear Flows (Ult) - Lateral Gust at 4000 
ft, pressure 9.4 psi, telescope brake off. 

4.4 

Beam Web Shear Flows - 2P condition, cabin pressure 
18.8 psi, no other loads. 

4.5 

Beam Web Shear Flows - Emergency Landing Condition 
9g forward. 

4.6 

Bulkhead Pressure Web Shear Flows - Lateral Gust at 
20,000 ft, telescope brake on. 

4.7 

Beam Cap Axial Loads, Forward Caps - Emergency 
Landing Condition 9g forward. 

4.8 

Beam Cap Axial Loads, Forward and Aft Caps, Half 
Bulkhead - Emergency Landing Condition 9g forward. 

4.9 

Beam Cap Axial Loads, Forward and Aft Caps, Half 
Bulkhead - Positive Gust at 20000 ft, brake 
engaged. 

4.10 

Beam Cap Axial Loads, Forward and Aft Caps, Half 
Bulkhead - 2P condition, 18.8 psi, no other loads. 

4.11 

Beam Cap Stresses, Forward Caps - IP condition, 9.4 
psi. 

4.12 

Beam Cap Stresses, Forward and Aft Caps, Half 
Bulkhead - IP, 9.4 psi, no other loads. 

4.13 

Pressure Web Shear Flows - Lateral Gust at 40000 
ft, brake off. 

4.14 

Beam Cap Areas Forward Bulkhead, Aft Face 

4.15 

Beam Cap Areas Forward Bulkhead, Fwd Face 
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4.0 STRUCTURAL DESCRIPTION AND ANALYSIS RESULTS (cont'd) 

4 . 2 AFT BULKHEAD 

The aft cavity bulkhead is located at FS 1960, and although of 
similar basic design, the structure is of less massive proportions. 
The four major functions of the bulkhead are; 

i) It supports the telescope brake and provides a reaction point 
for the brake axial and lateral forces. 

ii) It separates the telescope cavity from the unpressurized 
section of the aft fuselage. 

iii) It acts as the primary aft framing member for the fuselage 
cut-out . 

iv) It provides access, via a door, to the telescope cavity from 
the aft section of the fuselage for ground maintenance 
purpose . 

As stated above the basic construction of the bulkhead is similar 
to the forward bulkhead although in this case beams are only 10 
inches deep and are aligned vertically and horizontally only. The 
45 degree members are not required as the bulkhead does not have 
the 100 X 100 inch cut-out and the applied differential pressure is 
only 1.0 psi limit. A major consideration in the design of the aft 
bulkhead is stiffness, as it supports the telescope brake which is 
assumed to be a spring loaded device. The bulkhead may not deflect 
relative to the forward bulkhead to a point where the brake spring 
loses its effectiveness. Similarly, the bulkhead should not be so 
stiff that, when it is engaged and the forward fuselage 
pressurized, the brake creates excessive forces. 

The bulkhead contains a non-structural door. This door provides 
access to the cavity for ground maintenance and ensures a thermal 
seal during ground pre-cooling. 

The aft bulkhead pressure web is made from the same LS5150 T73 
material as that of the forward bulkhead with the risers facing 
forward into the cavity. The beams are also similar to those of 
the forward bulkhead, though only 10 inches deep and with caps of 
much less massive proportions. As with the forward bulkhead the 
beam webs are non-buckling up to limit load to obtain a constant 
spring rate for the telescope brake. This requires a 0.063 inch 
thick web with stiffeners at 5 to 6 inch spacing. 

As with the forward bulkhead, the beam caps are made from 7075-T73 
machined extrusion with similar joints where the vertical and 
horizontal beams cross. However, in this case due to the lighter 
loading and reduced criticality of this bulkhead, the caps are 
spliced using bathtub fittings only, with additional reinforcement 
at the door corners. 
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4.0 STRUCTURAL DESCRIPTION AND ANALYSIS RESULTS (cont'd) 

4.2 AFT BULKHEAD (cont'd) 

The recommended materials are the same as for the forward bulkhead. 

The results of the FEA model runs are summarized on the following 
pages as shown on the chart below: 


Page 

Description 

4.18 

Bulkhead Beam Bending Moments, 1.5 psi cavity 
pressure with Telescope Brake on. Section at WL 200 
and up between the Bulkhead Centerline and the Access 
Door Inboard Edge Framing Member. 

4.19 

Bulkhead Beam Bending Moments, 1.5 psi cavity 
pressure with Telescope Brake On. Section at WL 200 
and up from the Bulkhead Centerline outboard 
(opposite to Fig page 4.16). 

4.20 

Bulkhead Web Shear flows. Lateral Bending Envelope 
Condition. 

4.21 

Aft Bulkhead Properties. 
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4.0 STRUCTURAL DESCRIPTION AND ANALYSIS RESULTS (cont'd) 

4.3 BULKHEAD TO SHELL INTERFACE 

The interface between the forward bulkhead and the fuselage shell 
is a critical fatigue point on the modification for two basic 
reasons which are closely related: 

i) The fuselage is "tied across” due to the flat bulkhead and 
hence is constrained against expanding under cabin pressure 
loads as it does forward of the bulkhead. 

ii) The bulkhead being flat deflects in bending as opposed to 
stretching, as is the case with a membrane type bulkhead. 

These two phenomena tend to produce significant rotation between 
the bulkhead and the shell at the interface. To counteract this 
tendency, intercostals tailored to produce a gentle rate of change 
of slope, are used at the bulkhead beam to shell intersection. 
These intercostals are varied in depth and cross-sectional area, to 
produce relatively even rates of load change in the members and 
hence a good fatigue life. These intercostals were incorporated in 
the FEA model in gross terms except in the case of one intercostal, 
which for the purposes of this study was modified in more detail, 
and assumed to be typical. This intercostal was modeled to the 
point where individual fastener loads could be extracted and a 
Stress Severity Factor (SSF) calculated. This SSF was then 
factored up to account for the effect of surface finish, anodizing, 
and other deleterious factors to produce a Fatigue Quality Index 
(K) . This can be used in the same manner as a Stress Concentration 
Factor (Kt) to make a prediction of service life. As no spectrum 
of SOFIA usage is available the life estimate was determined using 
the stresses from cabin pressurization ratioed up by 1.15 and using 
this as a once per flight, (ground-air-ground) loading. This is a 
reasonable approach as pressurization loads are generally 
responsible for 90% of pressurized fuselage fatigue damage. 

The Stress Severity Factor (SSF) which is calculated using Lockheed 
Report SMN 389 is the ratio of the maximum stress in a joint to a 
reference stress and is a function of several factors as described 
below. A copy of an AIAA paper presented by Cornell and Darby 
describing the approach is given in Appendix A. 
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4.0 STRUCTURAL DESCRIPTION AND ANALYSIS RESULTS (cont'd) 

4.3 BULKHEAD TO SHELL INTERFACE (cont'd) 

Due to the nature of the joint under consideration the SSF equation 
reverts to the original equation developed by Jarfall at Boeing. 


SSF = Ktg X fbp = Ktb X 0 X fbr 

fref 

where = an empirical hole preparation and filling factor 

0 = fastener tilt factor 

Ktb = hole bearing stress concentration factor 
Ktg = gross area axial stress concentration factor 
fbr = uniform hole bearing stress 
fpb = gross area uniaxial hole by-pass stress 
fref = gross area uniaxial reference stress away from joint 

The critical point on the intercostal occurs at fastener 13 (ref 

page 4.24 and 4.25 for fastener locations and local loads.) 

Sample Fatigue Quality Index Calculation for fastener 13 

I A 



Ktg 

= 3.25 

Ktb 

= 1.26 

0 

= 1.09 


in 

• 

o 

II 


0.25 

A SECTION A A 

VIEW SHOWING HALF CAP AT FASTENER 13 
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4.0 STRUCTURAL DESCRIPTION AND ANALYSIS RESULTS (cont'd) 

4.3 BULKHEAD TO SHELL INTERFACE (cont'd) 13307 13325 
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4.0 STRUCTURAL DESCRIPTION AND ANALYSIS RESULTS (cont'd) 

4.3 BULKHEAD TO SHELL INTERFACE (cont'd) 


Elem Area Load Stress Fastener 

Load 


Beam 

23667 

0 . 49 

730.54 

1490.91 

94.89 


Beam 

23008 

0.49 

825.04 

1083.70 

94.50 


Bean 

23669 

0.49 

920.07 

1877.69 

95.03 


Beam 

23070 

0.4 9 

1017.99 

2077.53 

97.92 


Beam 

23671 

0.49 

1121.42 

2288.61 

103.43 

Fastener 

Beam 

23072 

0.49 

1228.04 

2500.20 

100.02 

Notation 

Beam 

23673 

0.49 

1337.43 

2729.45 

109.39 

Ref p 4. 24 

Beam 

23074 

0.49 

1450.20 

2971.90 

118.83 

1 

Beam 

23675 

0.49 

1582.25 

3229.08 

125.95 

2 

Beam 

23070 

0.4 9 

1207.85 

2405.0 0 

-374.40 

3 

Beam 

23677 

0 . 1 9 

881.69 

1799.37 

-326. J 6 

4 

Beam 

23078 

0.4 9 

519. (.0 

1000.41 

-302.09 


Beam 

23679 

0 . 49 

0.00 

0 . 0 0 

-519.60 


Beam 

23080 

0.7 5 

0 . 0 0 

0.00 

0 . f) (1 


Beam 

23681 

0.75 

1399.70 

1860.27 

1399.70 


Beam 

23082 

0.75 

2572.59 

3430.12 

1172.89 

5 

Beam 

23683 

0.75 

3712.41 

4949.88 

1139.82 

6 

Bca m 

23084 

(!. 75 

4998 . 00 

0004 . 08 

.1285.05 

7 

Beam 

23085 

1 . 00 

5-'.76 .65 

5470.05 

478.59 

8 

Beam 

23080 

1.00 

5900 . 81 

59(.0.81 

484 .10 


Beam 

23607 

1 . 00 

6 4 40. 01 

6446 . 01 

485.20 


Beam 

23088 

1.00 

0932 . 70 

0932.70 

480.75 


Beam 

2 3089 

1.00 

7419.05 

7419.65 

486 . 89 


Beam 

23090 

1 . 00 

7910 . 38 

7910 . 38 

490.73 


Beam 

23691 

1 . 00 

8420.40 

8420 . 40 

304 .08 


Beam 

23092 

1.00 

8921 . 05 

8921.05 

501.19 


Beam 

23693 

1 .25 

9436 . 40 

7549.12 

514.75 


Beam 

23094 

1 . 25 

9939 .14 

7951 . 31 

502 .74 


Bear. 

23695 

1 . 25 

10443 . 90 

8355 . J 2 

504.70 

9 

Beam 

23090 

1.25 

10970 . 00 

8781.28 

532.70 

10 

Beam 

23697 

1.25 

1 J 530.80 

9224.64 

554.20 

11 

Beam 

23098 

1 . 25 

12105.80 

9084.04 

575 . Of) 

12 

Beam 

23699 

1 .25 

12697.20 

lOJ 57.76 

59i . 40 

13 

Beam 

23700 

.1 .25 

12258 . 00 

9800.40 

-439.20 

14 

Beam 

23701 

1.25 

11375.70 

9100.56 

-882.30 

15 

Beam 

23702 

1.25 

10501.00 

8449.28 

-814 .10 

16 

Beam 

23703 

1.25 

9810 . 38 

7848.30 

-751.22 

17 

Beam 

23704 

1.25 

9111.87 

7289.50 

-098.51 


Beam 

23705 

1.25 

8455.67 

0764.54 

-650.20 


Beam 

23700 

1 . 25 

7831.44 

0205.15 

-024.23 


Beam 

23707 

1.25 

7236 . 25 

5789.00 

-595.19 


Beam 

23708 

1.25 

0058 . 04 

5320.43 

-578.21 
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4.0 STRUCTXJRAL DESCRIPTION AND ANALYSIS RESULTS (cont'd) 
4.3 BULKHEAD TO SHELL INTERFACE (cont'd) 


fref = 10157 psi 

fpb = 9864 psi 

fbr = 12442 psi 

SSf = 0.75 3.25 X 9684 + 1.26 X 1.09 X 12442 

10157 

= 3.59 

Fatigue Quality Index K = SSF x Surface Factors 

= 3.59 X 1.25 

= 4.48 

Increasing the reference stress by 1.15 to account for Flight Loads 

f = 11680 psi 

= 5840 + 5840 


From Lockheed SLM 16c Constant life diagrams 

K Life 

2 Infinity 

3 Infinity 

4 3 X 10E6 


6 9 X 10E4 

By interpolation from page 4.24 for K = 4.48 

Unfactored Life = 500000 GAG pressure cycles (flights) 

Factored Life = 167000 flights 

Transition stresses and shear flows for the 9.4 psi cabin pressure 
condition are shown on pages 4.28 thru 4.33. 

Pages 4.34 and 4.35 show the location of the transition 
intercostals and the bending moment diagram for the FS 1720 frame 
at the forward end of the intercostals. 




OJ o 



ESTIMATED FLIGHTS 

paws -062 



Input 
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4.0 STRUCTURAL DESCRIPTION AND ANALYSIS RESULTS (cont'd) 

4.4 CUTOUT SURROUND STRUCTURE 

The purpose of the cutout surround structure is to frame the cutout 
opening and restore the aft fuselage strength and stiffness lost 
due to the shell cutout, it does this by essentially acting as a 
Vierendel Truss using the primary and secondary frames and 
longerons to form I beams framing the cutout. This is shown on 
sheet 6 of drawing B45114000 and depicted schematically below. 



To permit operation of the cavity door the primary longerons (the 
ones immediately adjacent to the opening) are external to the shell 
mold line and faired with a fiberglass fairing for aerodynamic 
smoothness. The secondary longerons are internal. The skin 
between the primary and secondary members is doubled up to carry 
the high shear loads and a tripler is added to provide extra 
reinforcement to the skin immediately adjacent to the cutout. The 
longerons carry high axial loads and are unsupported columns 
between the bulkheads, a distance of 180 inches, and hence require 
large moments of inertia. 

Upper Primary Longeron required Moment of Inertia = 100 in4 
Upper Secondary Longeron required Moment of Inertia = 35 in4 
Lower Primary Longeron required Moment of Inertia = 70 in4 


Upper Secondary Longeron required Moment of Inertia 


50 in4 


All Longerons are assigned a cross-sectional area of 6 in2 for this 
study. 
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4.0 STRUCTURAL DESCRIPTION AND ANALYSIS RESULTS (cont'd) 

4.4 CUTOUT SURROUND STRUCTURE (cont'd) 

The longerons are of conventional construction with extruded or 
machined caps and a stiffened shear web between them. Where the 
caps must be spliced, to go through a bulkhead or frame, or to 
transition down to a reduced height, machined fittings are used. 

The stringers are spliced with fittings similar to the standard 
Boeing stringer splices where they pass through the two bulkheads. 
The basic difference is a flat web machined into the fitting. This 
web attaches to the bulkhead cap ring on the side opposite to the 
cavity. The flange prevents pressure leaks from the cabin into 
cavity where the fitting passes through the mousehole in the cap. 
A typical fitting is shown on sheet 8 of drawing B54114000. 

The cutout surround doubler at the forward bulkhead extends around 
the full circumference of the aircraft to form part of the bulkhead 
to shell interface described in section 4.3. The doubler also 
extends along the line of the longerons and cutout surround frames 
to aid in load runout. 

The tripler is of limited extent and made from segments of 2024 
plate, machine spliced together and tapered at the ends to aid in 
the longeron load runout redistribution. At the forward end of 
the cutout, the tripler extends from the cutout edge to the forward 
bulkhead cap. This provides added support and stability to the 
overhanging skin and door track frame. 

The preferred materials for the cutout surround structure are: 
Longeron caps 7075-T76 extrusion 

Stringer splices 7075-T73 hand forged billet 

Doubler 2024-T3 Clad 


Tripler 

Longeron Splice Fittings 


2024-T861 Plate 


7075-T73 bar 
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4.0 STRUCTURAL DESCRIPTION AND ANALYSIS RESULTS (cont'd) 

4.4 CUTOOT SURROUND STRUCTURE (cont'd) 

The results of the FEA runs are summarized on the following pages 
as shown on the chart below. 


Page 

Description 

4.43 

Upper Longerons Axial Load - Max Tension Load Scans 

4.44 

Upper Longerons Axial Load - Max Compression Load 
Scans 

4.45 

Lower Longerons Axial Load - Max Tension Load Scans 

4.46 

Lower Longerons Axial Load - Max Compression Load 
Scans 

4.47 

Stringer Axial Stresses around Cutout - Fuselage 
Downbendihg Envelope plus Cabin Pressure 

4.48 

Stringer Axial Stresses around Cutout - Fuselage 
Sidebending Envelope plus Cabin Pressure 

4.49 

Shell Axial Stresses around Cutout: 

Upper view Download Envelope plus Cabin Pressure 

Lower view Upload Envelope plus Cabin Pressure 

Left Side of Fuselage 

4.50 

Shell Axial Stresses around Cutout, Upper Fuselage: 
Upper View Side-load Envelope (right) plus Cabin 

Pressure 

Lower View Side-load Envelope (left) plus Cabin 

Pressure 

4.51 

Shell Shear Flows around Cutout - Side-load Envelope 
Condition 

4.52 

Shell Shear Stresses around Cutout - Side-load 
Envelope Condition 
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Stringer axial stresses around cutout 
Fuselage s>cie— bendi ng envelope plus cabin press-ur-e 
Re-f GSRR report 1391-1 page 21 
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4.0 STRUCTURAL DESCRIPTION AND ANALYSIS RESULTS (cont'd) 

4 . 4 FLOORING 

With the exception of a maintenance walkway opposite the cavity 
opening, the main deck floor is removed between the cavity 
bulkheads to permit installation of the telescope. The floor is 
not attached to either bulkhead except at the skin line to permit 
the bulkheads to deflect freely under the various load conditions. 
The cargo floor is retained as is except for being cut back from 
the bulkheads, and a maintenance walking surface is added. 

Forward of the pressure bulkhead and aft of FS 1600 the floor is 
cut away between BL SOL and 60R to permit installation of the 
telescope, and to provide comfortable access to the telescope 
instrument package. 

The shell of the 747 is not circular: when the shell is pressurized 
it wants to achieve a circular cross-section to most efficiently 
react the loading. This is normally resisted by the maindeck 
floorbeams which are partially removed. To simulate this restraint 
the cut section of the floor has an inner cap and shear web added 
at WL 199.3 to force it to act as a deep beam supporting the shell 
between the bulkhead and FS 1580. As with other sections of floor 
adjacent to the bulkhead, the floor/beam is only attached to the 
bulkhead at the skinline. 



Floor Schematic 
Ref Dwg B4511400 Sht 3 
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4.0 STRUCTURAL DESCRIPTION AND ANALYSIS RESULTS (cont'd) 

4 . 4 FLOORING ( cont ' d ) 

The results of the FEA runs are summarized on the following pages 
as shown on the chart below: 


Page 

Description 

4.55 

Main Deck Floor - Shear Flows at Added Shear Webs 
Condition - 2P (18.8 psi) 

4.56 

Main Deck Floor - Axial Forces in Members at Floor 
Cutout - Condition 2P (18.8 psi) 

4.57 

Main Deck Floor - Deflected shape of Floor Cutout 
Area - Condition 2p (18.8 psi) 

4.58 

Main Deck Floor Section Properties 

4.59 

Main Deck Floor - Added Shear Webs 

4.60 

Main Deck Floor - Bending Stresses and Displacements 
9.4 psi cabin pressure 
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ELEMENT T /PES 

RIGID 
: 5hF>^-{ 

Sl-0_L 



BERM: RBveie03 

GRID R: RIs=«C:5a63 
GRID B: RZW08ea4 
! Length= 30.0 
Xe: X= 1.3000 

r= .3000 
Z= .0003 
X= .030.3 
Y= -.'?SS0 
Z= .0536 
PBEHM: PDV 
fla= . ; 90 

Ila= . 3:3 
I2a= .0:3 
.310 



Zut-fliuay Cabin Floor Section 


Label 

Area 

11 

12 

J 



in' 2 

in'4 

in“4 

in"4 


1 

. 406 

. 10 

. 10 

. 10 

* 

2 

1.308 

7.619 

. 14 

.05 



1 . 155 

2 . 233 

.14 

.05 

* 

3 

2.996 41.499 

. 14 

.05 



2.656 60.421 

. 14 

.05 

* 

4 

1 . 315 

7 . 966 

. 14 

.05 



2.966 41.499 

. 14 

.05 


5 

1 .155 

2 . 233 

.14 

.05 



1.315 

7 . 966 

. 14 

.05 


6 

1.850 16.00 

3.00 

. 10 


7 

. 526 

. 10 

. 10 

. 10 

* 

8 

. 567 

.10 

. 10 

. 10 



. 526 

. 10 

. 10 

. 10 


9 

. 608 

.10 

. 10 

. 10 



. 567 

. 10 

. 10 

. 10 

* 

10 

. 649 

. 10 

. 10 

. 10 



. 608 

. 10 

. 10 

. 10 

♦ 

11 

. 661 

. 10 

. 10 

. 10 



. 649 

. 10 

. 10 

. 10 

* 

12 

. 673 

. 10 

. 10 

. 10 



. 661 

. 10 

. 10 

. 10 


13 

2.000 





14 

. 190 

.01 

.01 

.01 


* Tapered beam elements 
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ELEMENT T^'PES 

EZPII 

ROD 

FI-3ID 

ShFBn* 
9--ELL 


jriLE: cofr<^^ovJa^.s 

f 

LDCHSE: 2.0F?tESS 

XMIN= 1540 
XMRX= 184C 

YMIN= e 
YMHX= 3ee 

ZMIN= 199 
ZMRX= 20Q 


=LF"If 5 — Floor Fortuard of Fuel Bulkhead 
Elsforjned Shape Magnified 50 x 
E’erftent Types and Connectivity 
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f2iaZO<32C 
RI80OG2 & 
I^I'7>SoQ2€- 
R.I7&0G2€' 
RI7<^0GZ€ 

R.nzoG2e, 
fil700GZ€ 
R./6SoGZ€i 


Displacements 


Grid 

ID 

X 

(in) 

Y 

(In) 

z 

( in) 

R1660G26 

0.767495 

0.098383 

-0.885930 

R1680G26 

0.774250 

0.101986 

-0.908760 

R3700G26 

0.780700 

0.084695 

-0.925645 

R1720G26 

0.786195 

0.043044 

-0.934510 

R1740G26 

0.791385 

-0.039607 

-0.932870 

R1760G26 

0,795875 

-0.047178 

-0.951515 

R1780G26 

0.797060 

0.060308 

-0.994550 

R1800G26 

0.795190 

0.058998 

-1 .016920 

R1820G26 

0.794140 

0.046490 

-1 .032235 


/ZIC 60 G^£i 


Grid 

Ml 

( in-lbs) 

12 

(ln‘4) 

c 

( in) 

Bending 

Stress 

(psi) 

R1660G26 

2639.0 

16.000 

— 

— 

R1680G26 

11371.0 

16.000 

— 

— 

R1700G26 

11264.0 

16.000 

— 

— 

R1720G26 

12792.0 

2.233 

1 .389 

7957.05 

R1740G26 

15458.0 

7.966 

2.662 

5165.60 

R1760G26 

180766.0 

41.499 

3.935 

17140.51 

R1780G26 

59586.0 

7.619 

2.604 

20365. 13 

R1800G26 

2416.0 

2.233 

1.389 

1502.83 

R1820G26 

253.0 

0.100 

— 

— — 
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4.0 STRUCTURAL DESCRIPTION AND ANALYSIS RESULTS (cont'd) 

4 . 6 STIFFNESS 

4.6.1 Overall Airframe Stiffness 

The relative stiffness of the baseline 747-100 and the 747 in the 
SOFIA configuration were compared by applying unit loads to the aft 
fuselage at FS 2180 with the fuselage being held fixed at FS 1480. 
The loads were applied to the shell through rigid elements which 
force the shell to hold its shape at the loaded end. 

The results of this stiffness comparison are shown below: 


Unit Load Type 

Lateral Py 
Vertical Pz 
Torsion Mx 

4.6.2 Bulkhead Stiffnesses 


SOFA Stiffness/ 
Baseline Stiffness 

1.025 

1.019 

0.903 


The deflection of the forward (pressure) bulkhead under 9.4 psi 
cabin pressure is shown on page 4.51. The deflections shown are 
relative to an unloaded airframe. The relevant deflections are 
those of the telescope attach points at the four corners to the 
shell. This deflection is 0.531 inch maximum at the lower right 
hand corner. 


The deflection of the aft bulkhead is shown on page 4.52. The 
condition shown is for 9.4 psi cabin pressure and 5000 lb. limit 
brake load. The cabin pressure has no effect on the bulkhead 
directly but does serve to move the shell datum aft for comparison 
with the forward bulkhead deflection. No conditions were run with 
the bulkheads coupled through the brake spring as its 
characteristics are currently undefined. 
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4.0 STRUCTURAL DESCRIPTION AND ANALYSIS RESULTS (cont'd) 
4.6 STIFFNESS (cont'd) 

4.6.2 Bulkhead Stiffnesses (cont'd) 



Forward (Pressure) Bulkhead Deflections 
9.4 psi Cabin Pressure 
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4.0 STRUCTURAL DESCRIPTION AND ANALYSIS RESULTS (cont'd) 
4.6 STIFFNESS (cont'd) 

4.6.2 Bulkhead Stiffnesses (cont'd) 


d.'527 





Aft Bulkhead Deflections 
Condition 9.4 psi plus Brake Load 
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5.0 WEIGHT AND BALANCE 

5.1 WEIGHT ESTIMATE 

The weight of the SOFIA structure covered by this report. Forward 
Bulkhead, Aft Bulkhead, Bulkhead to Shell Interface, Cutout 
Surround Structure, and Flooring has been estimated using the FEA 
model with the weight change calculated for the before-mod and 
after-mod conditions. The structural weight estimate makes no 
allowance for Telescope local supports which are part of the 
Telescope/Bearing assembly nor does it reflect local structure 
required to. support the Cavity Door, its associated mechanisms and 
equipment support. 

Weight change of aluminum structure (from FEA model) 

= + 7500 lb 

Adding 35% for fasteners 
and stiffeners, etc. 

15% contingency 


+ 10125 lb 
+ 11640 lb 
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6.0 FUTURE TASKS 

The cavity structure, bulkheads, surround structure, and bulkhead 
to shell transition should be investigated more deeply to validate 
the design concepts expounded here. This investigation should be 
down to the detail stress analysis level to ensure that design 
items accepted as feasible under the limited budget of this study 
can, in fact, be implemented. 
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CORRELATION OP ANALYSIS AND lEST DATA 
ID IHE EFFECT OP PASTEI^ER LOAD TRANSFER ON FATIGUE • 

B. L. Cornell and L. G. Darby* 

Lockheed-Georgia Company 
Marietta, Georgia 


Abstract 


Presented is a Stress Severity Factor method 
of fatigue cuialysis for mechanically fastened, 
single shear alunrinum lap Joints, with one or two 
rows of titanium fasteners. Also presented is 
correlation of analysis predictions to test data 
from two full scale wing fatigue test articles 
representing the C-3A trans 3 )ort. Fatigue damage is 
calculated for each loading condition recognizing 
that the stress severity factor varies within each 
load cycle; considering parameters of fastener load, 
uni ATd stress, diameter to thickness ratio, fas- 
tener head configuration; and using joint 5-N data 
which reflects fretting. Damage is simned for a 
spectrum using Palngren-Miners taeory. 

1 ,0 Introductian 

The significance of load transfer in mechan- 
ically fastened structural Joints naa oeen increased 
in recent years by the development of high strength 
fastener systems in the ils; TTS range, 

This development has tended to aaae structural 
Joints more efficient statically, out not more ef- 
ficient in fatigue. The desi^. oileana of creating, 
a balanced structural Joint, wnicr. satisfies both 
the static and fatigue criteria, is tnt establish- 
ment of the design hole beanr.g stresses. With 
strict design weight controls and tr.t availability 
of high strength fasteners, it is natural for the 
designer to want to exploit the full fastener 
strength capability, thereby intrccucing high hole 
bearing stresses which may have an adverse effect 
in fatigue. To provide a technical oasis to deter- 
mine the design hole bearing stresses requires a 
fatigue analysis method which recogrures the 
effects of fastener load transfer. With such a 
method the structural Joint can be designed for 
fatigue, and with the well established static 
methods of analysis, a balanced joint design can 
be obtained. 

This paper will address the effects of fastener 
load transfer on the fatigue behavior of a single 
lap shear Joint. A fatigue analysis method is 
presented and the results compared to existing 
methods of analysis using the test data obtained 
from two full scale fatigue test articles. 

2,0 Fatigue Analysis Method 

In a practical fatigue analysis method, for 
any type of mechanically fastener Joint, it is 
advantageous to use existing technology as much as 
possible. The fatigue analysis method presented 
for single ATngrinmn lap shear Joints is, therefore, 
based on the conventional i^-method of a n a l ysis. 

Ihe readily available empirical constant amplitude 
notched S-N data are utilized, but are modified to 
Joint 5-N data by accounting for hole fretting 
effects. Fbr Joint damage calculations the well- 
known PalmgrenAtlner’s linear cumlative damage 

* Aircraft Development 
Fhgineer Specialists 


theory is used. 

Damage Equation 

The Palmgren/^liner ’ 3 equation can be expressed 
as follows: 



Where 

Dm ■ Total Damage, failure occurring when 

« Applied cycles at stress level m 

K « Allowable cycles at stress level m 
m 

In this equation the value N is obtained 
from the S-N data and is a ihnction of the uniaxial 
stress cycle and the associated Joint stress con- 
centration factor (Km)« Since the uniaxial stress- 
es and number of applied cycles (n ) are always 
Imown, the solution for the total damage depends 
on the determination of the Km values. 

In the past the Krp present in a Joint has, in 
general, been determined experimentally. Jarfall^M 
however, investigated the problem of analytically 
defining the significant parameters affecting the 
stress concentration factor in a lap axial. Joint, 
and proposed the "Stress Severity Factor (SSP) 
Method". Since SSP can be equated to K-f the SSP 
Method is used to determine the Joint Km* The K-- 
Method of fatigue analysis is, therefore, redefined 
and called the SSP-Method of fatigue analysis. 

Stress Severity Factor Equations (SSP) 


The basic equation for the SSF proposed by 
Jarfall for a lap axial Joint is as follows: 

(UipISal)' 

Where • Gross area uniaxial reference 

stress away from Joint (psi) 

f B Gross area uniaxial hole by-pass 
stress (psi) 

f. • Unifonn hole bearing stress (psi) 


Sg 


B Gross area uniaxial stress concen- 
tration factor 

B Hole bearing stress concentration 
factor 


6 B Fastener tilt factor, ratio of uni- 
fonn to peak hole bearing stress 


a 0 m ^nplrleal hole preparation and hole 
filling factors, respectively 

Unfortunately.Jthis SSP equation, in its pres- 
ent form, cannot be used for a lap shear Joint as 
the load path of the transfer load is different. 
Figure 1-0 compares the load transfer characteris- 
tic of the lap axial and lap shear Joints. It can 
be seen that, in the lap axial Joint the load 
transfer produces a uniajo.al stress gradient across 
the joint with no shear stress gradient, while in 
the lap shear joint the exact opposite is true. 

The load transfer in the lap shear Joint produces 
a shear stress gradient across the Joint and no 
uniaxial stress gradient. 


To account for these affects, and make the SSF 
equation applicable to a lap shear Joint, tiie basic 
S3F equation is modified expanded into the fnnn: 


SSP . 

(Lap Shear) ^ax 


^ ® "br 



Where ■ Gross area uniaxial stress acting 

along the Joint (psi) 

T » Gross area shear stress away from 
Joint (psi) 

» Gross area shear stress concentration 
factor compatible with 

The remaining terms being as defined in the SSF 
equation for the lap axial Joint. 


in both the lap axial and lap shear SSF equa- 
tions, the terms inside the brackets describe the 
magnitude of tne peak stress in the Joint material 
at the edge of the fastener hole. Figure Z,j shows 
a pictorial breakdown of the peak stress into the 
three stress components present in tne lap shear 
Joint. Ihe first and third terms are, simply, 
gross area stresses times a stress concentration 
factor and give tne peak stress components due to 
tne tuiiaxial and shear stresses, respectively. 

The second, or middle, term is similar, and gives 
the peak stress con^xjnent due to the hole bearing 
stress. However, the term contains the fastener 
tilt factor (B), By definition, B is equal to the 
ratio of uniform to peak hole bearing stress, and 
is req^ured to account for the peaking that occurs 
ax the faying surface of the Joint when the fasten- 
er deforms (tilts) under the action of the load 
transfer. The value of B depends on whether the 
Joint is of the single or double shear type. An 
approximate analytical solution for 0 , in a single 
aluminum lap shear Joint, is presented in Section 
4.C. Curves of ^ versus X./D for two thictaiess 
ratios using titanium fasteners are also presented 
in this section. To derive values for the stress 
concentration factors K^, K^, and a typical 
single lap shear Joint vath two rows oT fasteners 
is analyzed using finite element techniques. A 
discussion of this analysis and final values for 
the stress concentration factors are presented in 
Section 3«0. 

To account for hole preparation and hole 
filling effects the peak stress is corrected to an 
apparent peak stress by applying the a ^uld 0 
factors. These factors are empirical and must be 
determined by test. The finite element analysis 
of Section uses the condition of a perfectly 
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finished hole and a neat fitting fastener as the 
baseline condition. A hole and fastener instal- 
lation which meet these conditions could, 
toerefore, be expected to demonstrate a = 1 and 
3 » 1 . For holes that have a less than perfect 
finish the a factor should increase and for fas- 
teners with interference or clamp-up the 0 factor 
should decrease. The approach used to detennine 
the a and 0 factors on two full scale test 
articles is discussed in Section 5.0. 

The final step in determining the SSP is to 
divide the product of a , ^ and the peak stress by 
the gross area uniaxial stress (f ). 5y perform- 
ing this operation it is apparent^^hat, if f 
approached zero and the terms containing f. ^and 
T had some finite value, SSF would approach infin- 
ity. In a structure, such as a wing box, this con- 
dition could only be achieved if the wing is 
subjected to pure torsion. In a realistic loads 
environment a wing will never experience Just pure 
torsion; some bending must always be present. Hence 
this extreme condition should never arise. 


LAP AXIAL JOINT LAP SHEAR JOINT 



FIGURE 1.0 LOAD TRANSFER - LAP AXIAL AND 

lap shear joints 


GROSS AREA UNIAXIAL STRESS 



FIGURE 2.0 PEAK STRESS COMPONENTS - 
LAP SHEAR JOINT 


Uniform Hole Bearing Stress ^nd Shear 

Stress ( T ) 

For a single lap shear Joint, it is convenient 
to express the fatigue loads in terms of the gross 
area xiniaxial stress and shear flow (q). Again, 
referring to Figure 1..^ which shows the load trans- 
fer characteristic of the lap shear Joint, the 
fastener loads can be readily calculated using the 
equation: 


4nd therefore 


^br - “ <>/-. 

'Where Ap • Pas^eIle^ Load (lbs) 

a » Fastener pitch (ins) 

a » Number of rows, i or 2 rows 
only 

D « Fastener Diameter (ins) 

t > Thicioiess of Joint element 
to be analyzed (ins) 

For the case n « 2 and unequal thlclmess 
elements, the fastener loads become statically 
indeterminate. However, providing the thictoess 
ratio is not less than 0.7, analysis shows, the 
error in using the above equations is small. 

Allowable Cycles (N) 

Knowing the applied stress cycle and the SSF, 
it is a simple matter to enter the- S-N data and 
determine the allowable cycles (N). However, in 
a lap shear joint where the shear flow induced 
load transfer is independent of the gross area 
uniaxial stress, the SSF will vary throughout the 
loading cycle. Also, if the gross area uniaxial 
stress cycle is tension-cos^ression-tezxsion, the 
SSF will vary due to the different stress concen- 
tration factors for tension and conqjression. For 
example, taking the later ease of a tension- 
compression- tension gross area uniaxial stress 
cycle with no load transfer, the stress concentra- 
tion factors (K-^), from Section 3*0 for hole (1), 
are 3-2 for tension and ^ .35 f'or compression. If 
a and 0 are equal to t.O, SSF « 3-2 and 
SSP - 1.35. 


To account for a variable SSF, a numerical 
solution is developed based on the observation 
that within a given S-N data set it is possible to 
find a single equivalent value of SSF, which, when 
used with the gross area uniaxial stress cycle and 
applied to the apparent peak stress cycle, the 
allowable cycles (N) remain a constant. Intro- 
ducing the S-N data mean (m) and variable (v) 
stress notation, the numerical solution may be 
expressed as: 


N for 





N for 





Where SSF^ » Equivalent SSF 


As part of this solution it is necessary to 
rearrange the lap shear SSF equation Into the form 








S f, 


br 




3SP 


W, 


az 


Where f 


pk 


apparent peak stress 
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A sinqple step-by-step procedure to the solu- 
tion for K and SSF^ is as foliov;s; 

Given the fatigue condition 


Snax 


“ax . ’ %iin 
min . 


n cycles 


1 ) Determine f and f 

ax,^ ax 

V a 

2) At the desired analysis location within tne 
joint determine 


pk 

m 

3) Detemine f 


and 




“pk . 
min 

and f 

pk. 


m 


4) Determine SSF. « f . and SSF« - f ^ 

f f 

ax„ ax 


5) If SSE.^ is equal to SSF 2 , enter S-N data with 

^ax ' ^ax ^ aetermine allowable 

cycles (N)- 

SSF- « SSSP. » SSP^ 

6) If SSF^ is not equal to SSF- enter S-N data 
with f^ , f^ and SSF. to determine allow- 

able cycles N . Ehter S-N data with f . /SSF. 


1 

f ^ /SSP^ and SSF^ to determine 
allowable cycles 

If > N 2 increase SSP^ until ■ 
If < N 2 decrease SSF^ until 
Allowable cycles (N) a N. « 






SSF^ 


SSF^ at 




Figure 3-0 presents typical curves of Nj and 
N 2 versus SSF^ and clearly shows that there is 
only one unique solution for N and SSF-. 



Under complex dynamic loading it is possible 
for the f and q cycles xo be out of phase. In 
such a cate, the solution should be repeated for 
different time intervals to find the most critical 
values for N* and SSP^. 


material 7075- T6 PLATE 




< ALLOWABLE CYCLES 

^ FIGURE 4.0 effects OF FREHInG ON notched 
SPECIMEN FATIGUE auOw*8LES 


Data 


5 -N data suitable for fatigue analysis of 
mechanically fastened joints oust be based on 
joint data. Also, to be compatible with the S3F- 
method of fatigue analysis the data oust be 
expressed in terms of SSF or K». 

In Figure ^.C the results of Fatigue Tests 
on notched specimens, with and without fretting, 
are presented, ( 2 ) The specimens are aacmned from 
7075- T6 plate ma^rial and notcned to a K„, « 5-25* 
The figure shows the reduction in. fatigue ^allow- 
ables in the low stress-high cycle region of the 
S-N curve, when fretting is present, is consider- 
able. At 10^ cycles, for example, the fatigue 
allowables are reduced by ^ 0 $. 

Since fretting does significantly change the 
shape of the S-N curve and cannot be avoided in 
joints, the trends of Figure 4,0 are used to modify 
the readily available notched S-N data for 7073- T6 
plate material. Ihe modified notched S-N data are 
considered applicable to joints and are called 
joint S-N data. By taking this approach the re- 


MATL 7075- T6 PLAH 
mean stress = 10 KSI 



FIGURE 5,0 JOINT S-N DATA BASED ON NOTCHED 
DATA WHEN FRETTING 
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suiting joint S-N data are still expressed in terms 
of SSF or K.„ and can be used with the SSP-method 
of fatigue analysis. 

To provide a consistent method to modify the 
notched S-K data to joint S-N data it is argued 
that the cycle relationships' of Figure 4.C between 
the notched curve and the notched curve with fret- 
ting, at any given value of allowable cycles on 
the notched curve, will remain a constant for all 
values of K_. Therefore, the two curves of Figure 
4.0 can be rcplotted in the form verstis 

Kjjotched* equations to fit this curve 

will give: 


Branch 

( 1 ) 

LOG 

N, . = 

Joint 

LOG K ^ ^ 

notched 


Branch 

( 2 ) 

LOG 

N . . 

joint 

.375 LOG N ^ ^ ^ 
notched 

+ 3.0 

Branch 

(3) 

LOG 

N. . 
joint 

.196 LOG i) ^ ^ 

notched 

+ 4.03 


These equations are used to compute values of 
joint ^"joint approximately 10 ^ cycles. 

At this point the qurves are faired to intercept 
the abscissa at 10^ cycles. Ebcamples of the joint 
S-K data generated in this manner for 7?75-T6 
plate material are presented in Figure 5.0. 






UNIT 

BENDING 



FIGURE 6.0 MODEL A - SINGLE LAP SHEAR JOINT 


3.0 Derivation of Stress Concentration Factors 


To derive the stress conoentration factors, a 
typical spar web to cap single lap shear joint with 
two rows of fasteners is analyzed using a con^t- 
erlzed linear finite element displacement method 
of analysis. The structural features of this 
typical Joint are represented in two mathematical 


models. Model 'A* azid Model 'B' « A brief descrip- 
tion of these models and their method of loading 
is given below. 

Finite Element Models 

The stzxicture considered in Model 'A* is shown 
in Figure 6.0; in this model the web, cap and 
stiffeners are simulated, using a coarse grid net- 
work of isotropic triangular membrane plate ele- 
ments and the fasteners are simulated using linear 
shear springs. Two sets of external loads are 
applied, a unit bending coxidltion with an uniaxial 
stress of 10 ksi at the cap riser decaying to 5,5 
ksi at the edge of the web* and a unit shear con- 
dition with a constant shear flow cf 1 OOO Ib/in 
in the web. The intent of Model 'A* is to provide 
boundary loads (inplane stresses and fastener loads) 
for the more detailed Model *B' . 

The riser detail, shown hatched in Figure 6,0 
and considered In Model *B* , Is presented in Figure 
7.0. The riser detail is sisulated using the same 
type of element as used in Msdel 'A' , except in 
M 5 del ’B' the grid size is reduced. The grid size 
around the two holes of interest, holes (I) and (2), 
in particular, is reduced to an extremely fine grid 
for mAT-iTw^m accuracy. In Model '3' the influence 
of a varying fastener bearing area with load is 
also recognized. A special program is used which 
considers the fastener as a series of rigid radial 
spokes, each spoke being connected to a grid point 
around the periphery of the hole. Tne fastener 
load is applied at the hub of these radial spokes. 
Since the actual fastener cannot transmit radial 
tension loads, the radial spokes must also satisfy 
this condition. “Ihe program is therefore arranged 
to iterate until all tension spokes are eliminated 
and only compression spokes remain to react the 
applied loading. It should be noted that a neat 
fitting fastener is the reference condition and no 
fastener interference is considered in this analysis 


3/16 OlA FASTENER -NEAT FIT (TYPICAL) 
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SIGN CONVENTION - POSITIVE LOADING SHOWN 

NEGATIVE LOADING OPPOSITE 

FIGURE 7.0 model B RISER DETAIL 

To load Model 'B* the unit boundary loads from 
Model 'A* are factored and added to form combined 
cases. Twenty-sever (27) combined load cases are 
considered and cover a range of uniaxial riser 
stresses from -20.0 to 3C.0 ksi with shear flows 
from 0 to 3728 Ibs/in. To be consistent with the 
actual wing combined loading, the shear flow is 
taken negative when the uniaxial riser stress is 
negative. The sign convention used throughout the 
analysis is shown in Figure 7*0. 

Finite Element Analysis Results 

In interpreting the results of the finite 
element analysis, tangential stresses are used 


COMBINED TENSION AND SHEAR 

HOLE M) HOLE ( 7 ) 
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Z 

Uj 

o 

z 

< 


f M = 0.5 . L / f = 3.0 

br ox br ox 


COMBINED COMPRESSION AND SHEAR 



. . . i . » . 1 1 -1 1 i 1 ^ 1 * 

90 120 240 270 300 60 90 120 240 270 300 

<6 (DEGREES) 

FIGURE 8.0 tangential STRESS DISTRIBUTION, LOADED HOLES 
IN UNIAXIAL STRESS FIELD OF 10 KSI 


InsteAxl of the more coaoon nwxlnigri and mnlmm 
principal stresses. Hiis approach is taken because, 
in fatigue analysis, it is a requirement that the 
stresses acting on-^onoon planes be con^ared and 
used to con^te damage, the use of tangential 
stresses ensures this requirement Is met. Further- 
more if it is assumed the tension loading is the 
most important loading component in fatigue, no 
significant error is introduced as the' tangential 
stress and principal stress are, in essence, 

the same in the highly stressed region of the hole. 

Comparison of the tangential stress distribu- 
tions with increasing shear flow clearly show the 
significance of fastener load transfer. Figure 3.0 
illustrates this trend and shows the tangential 
stress for holes (i) and (2) as a function of angu- 
lar location ^ . The values plotted are for a 
uniaxial tension and compression riser stress of 
1C ksi with two values of fastener load transfer, 
f /f equal to 0.5 and 3.0. It should be noted 
tglt.Ho relate the finite element results to the 
fatigue analysis method, the uniaxial riser stress 
is now taken equal to the gross area uniaxial 
stress, f^, defined in Section 2.0. 

Prom Figure 8.C it is iniDediately apparent the 
behavior of the tangential stress in compression is 
considerably different to the more classic behavior 
of the tangential stress in tension. Fbr example, 
at hole ( 1 ) , edge A, the peak values in compression 
are in the range C.5 to 1.2 times the uniaxial 
stress; while the peak values in tension are in the 
range of 3*6 to 5.9 times the uniaxial stress. 

Also, the figure shows the angular location of 
these peak values do not coincide, the angular lo- 
cations for both tension and comprtasxon being a 
function of the amount of fastener load transfer or 
shear flow. 

A fatigue analysis method, which recognizes 
the tangential stress distribution as a variable, 
is extremely cumbersome. Therefore, in the deriva- 
tion of the stress concentration factors the - 
following simplifying assumptions are made: 

( 1 ) Under combined tension and shear loading it 
is assumed the peak tangential stress is 
equal to the calculated peak tangential stress 
and no correction is applied for location. 

This is conservative, but not overly so, as 
the arc formed by the points of the calculated 
peak tangential stresses, and the stress 
gradients within this arc, are small for the 
range of loading investigated. Figure 9.0 
shows the variation in the angular location 
of the calculated peak tangential stress as a 
function of the ratio. The arc shown 

is approximately. So* for the range 0 

to 4.0. 

(2) Under combined compression and shear loading 
it is a requirement for fatigue analysis that 
the tangential stresses act at the same loca- 
tion as the peak tangential stresses for 
combined tension and shear loading. It is 
assumed that this requirement can be adequately 
satisfied by using average values of tangen- 
tial stress actiiig within a conanon arc. In 
the analysis, the arc formed by the points of 
the calculated peak tangential stresses under 
combined tension and shear loading is used. 
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0 1.0 2.0 3.0 4.0 

f ( 
br Ox 

FIGURE 9.0 VARIATION IN ANGULAR LOCATION OF 
PEAK tangential stress 



FIGURE 10.0 variation OF PEAK STRESS RATIO 
COMBINED TENSION AND SHEAR 


Applying a5Suinptior.s (1) and (2) to the finite 
element results yield the curves presented in Fig- 
ures and 11. 1. In both figures the variation 

of the peaK stress raxio, f , is piotoed 

versus the ratio, Snerl^f is the peak 

tangential stress. It trlould be n8€ed that, in 
Figure 11.0, f . is not the ziue value of the peak 
tangential stress, but the average value '/h.thin 
the considered comnon arc. Also, as the tangential 
stress under uniform bearing is always positive 
(tension), irrespective of the direction of the 
bearing load, f^^ is always taken positive. 

Combined Tension and Shear Case 

Figure li;.v gives the variation in the peak 
stress ratio for the combined tension and shear 
case. The results are linear, but show signifi- 
cant differences in slopes, the slopes increasing 
when moving from hole (1) edge A to hole (2) 

Edge a. This difference in slopes is the influence 
of shear stress on the peak tangential stress and 
is the reason for the expansion in the stress 
severity equation discussed in Section l.C. The 
curve for Hole (2) Edge A is omitted, as this 
location is not critical. 

The equations which describe these cxirves, 
rearranged to give the peak tangential stress, 
are as follows: 
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Kole (1) Edge A. = 5.2 , 0.9 

Hole (1) Edge 5, = 3.2 * 1.3 

Hole U) Sdge-B, = 2.i . 2.4 

Now, using the arguments that the shear 
stress component of the peak tangential stress at 
hole (1) Edge A is small, and the beazhng stress 
component is constant at all locations, it is 
possible to separate the beax*ing and shear stress 
components of the peak tangential stress. Pro- 
ceeding in this manner, the equations may be 
rewritten in the form; 

sole (1) Edge A. = 3-2 ♦ 0.9 

Hole (1) Edge 3, - 3.2 0.9 0.4 

Hole (2.) Edge B. - 2.6 0.9 1-5 f^r 


Ihe additional terms being the snear stress 
component expressed as a function of the bearing 
stress. Prom the finite element results the 
bearing stress can be redefined in terms of shear 
stress (T), Making this substitution and taking 
Edge 3 as the critical location, tn# final equa- 
tions for holes (l) and (2) for tne combined 
tension and shear case become: 


Hole (1), fp^ « 3.2 ♦ 0.9 f^j. 


Hole (2), f . o 2.8 f 
' ^ ’ pit ax 


*br 


0.64 T 
2.2 T 


Introducing the stress concentration factor nota- 
tion of Section 2.0 


Hole (1), . 3.2, - C.9. - 0.84 

Hole (2), . 2.8, . C.9. - 2.2 

Combined Compression and Shear Case 

Figure 11.0 gives the variation in the peak 
stress ratio for the combined coopraasion and shear 
case. Unlike the curves for combined tension and 
shear, these curves are non-linear, but exhibit the 
same differences in slopes due to shear stress. 


Also, with increasing values, the 

curves cross the abscissa showinf^that the tan- 
gential stresses, within the arc of interest, 
reverse and go into tension if the ratio of 


Ihe equations which describe these curves, 
rearranged to give the peak tangential stress, are 
as follows. Two sets of equations are required: 

Range -6.0 $ 0-'^ 

Hole ( 1 ) Edge A, 

V “ '-°53 - (-.055 .087) 

Hole (1) Edge B. 

V ■ .263) 


Hble (2) Edge 5, 


PK 


1.35 f - (-.2'.'7 , > .763) f, 


br 


Hange - « < 


-6.J 


Hole 

(1) 

Edge 

A, f . s 

pic 

4.0 

f 

ax 

:.9 

‘br 

Hole 

(1) 

Edge 

pk 

4.: 

f 

ax 


"br 

Hole 

(2) 

Edge 

B, f ^ = 

pk 

4.0 

f ♦ 

ax 

2.4 

"br 


Applying the same arguments as used for the 
combined tension and shear case, the final eqxia- 
tions for Holes (1) and (2) for combined compres- 
sion and shear become: 


Range -6.0 g S ^-0 

•Hole (1) 

'pk = ^-55 * (-.055 * .087) 

» (-.034 * .363) T 

Hole (2) 

V = * (-.055 * .087) V 

« (-.221 * .934) T 

Range - oc < < -6.0 

Hole (1). f . . 4.0 f 0.9 T -k 0.64T 
pk ax br 

Hole (2), fpj^ - 4.0 . 0.9 . 2.2 r 

Again, introducing the stress concentration factor 
notation of Section 2.0 

Range -6.0 g 0.0 — < < -6.0 

Hole (1) 

Kro 1-35 4.0 

(-.055 * .087 ) 0.9 

K,g (-.084 ♦ .368) 0.B4 

Hole (2) 

1.35 4.3 

Kj3 (-.055 * .087) C.9 

Krs (--22' W^ax * -5®^) 2.2 



FIGURE 1 1 .0 VARIATION OF PEAK STRESS RATIO 

COMBINED COMPRESSION AND SHEAR 
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Application 

To apply the stress concentration factors to 
a single lap shear Joint requires an \mderstanding 
of the Joint geometry and internal load distribu- 
tion. If for instwce, holes (1) and (2) shown 
in Figure 7.0 are considered to be in rows (1) and 
(2), respectively, the riser detail is considered 
as element ( 1 ) and the attacning web is considered 
as element (2). The stress concentration factors 
presented can be subscripted as follows: 


Element (l). Row (l); ^Sl i 

Element (l), Row (2); Sbi2’ ^12 

El em e n t (2), Row (l); K,jq2i ’ ^4321’ ^4s2l 

Element (2), Row (2); KrpQ22* ^4 b22’ ^4s22 

Now, to maintain the same geometry and shear 
stress distribution as element (i), the hole loca- 
tion in element (2) must be interchanged. Hence: 


S?G22 ■ ^11* ^4b 22 * ^rS22 “ ^S11 


^21 * Sgi2’ ^4b21 " ^12’ '4s2i “ ^TS12 

Fbr a lap shear Joint witn one row of fasten- 
ers the stress concentration factors with sub- 
scripts (11) are used for elements (1) and (2). 
However, to account for the iugner snear stresses 
a factor of 2.0 is applied to 
Limitations 


The finite element results presented are 
based on the following geometneal ratios. 

e/T) « 2.2, a/D - 4,9,. b/D • 3-6 

Where e is the edge distance, a is the pitch, b 
is the distance between rows and D tnt fastener 
diameter. However, further analysis %/ith the 
fasteners increased to l/4 inches in diameter, 
show no significant chaxige in the results. The 
lower limit for the geometrical ratios investigated 
can therefore be taken as: 

e/D * 1.68, a/t) - 3.7, b/D - 2.76 

Using the above ratios it is considered the 
stress concentration factors presented can be 
safely applied to single lap shear Joints which 
fall within the range 


1.5 

% 

e/D 


2.5 

3.5 


a/D 


5.0 

2.75 


b/D 


4.0 


4,0 Analytical Fastener Tilt Factor (^) 

As previously stated in Section 2.0, the ac- 
tion of the load transfer causes the fastener 
system to deform (tilt) and peaking of the bearing 
stress- occurs at the faying surface of the Joint. 
This peaking of the bearing stress is a function of 
the Joint configuration and will vaiy depending on 
the material of the plate and fastener. Joint 
thielmeaa to fastener diameter ratio and type of 
fastener head. 


To derive a linear solution for the peak 
bearing stress at the faying surface and, therefore, 
the tilt factor (^) to be applied to the unifer® 
bearing stress, requires an extremely sophisticated 
method of analysis. In addition to the usual plate 
stifihess, fastener bending and shear stiffness 
tezms, the method should include tezms to account 
for the tension load introduced into the fastener 
system by the clamping action of the fastener head. 
HazTls, OJalvo and Hooson^^) proposed a linear 
solution for the peak bearizxg stress, in a single 
lap shear Joint, cozisidering the fastener as a 
shoi^ beam on an elastic foundation. The solution, 
however, did not include fastener tension effects 
and is^sed the arbitrary boundazy conditions of a 
negligible-head or rigid-head. To expand their 
solution to include the tension load effects would 
be a formidable task and is not attempted. Instead, 
to derive the equatiozis for the tilt factor. 



a P . M 


p(3) H 




I 



AP.Mp^ 


1*^ bf p® 

br(3D 


f 

p® 

PLATE BEARING FORCES 


^ 



I 

! 




f I 


FASTENER TENSION FORCES 


M =AP/^ (L M /.V * ^ mk 

'2' \ r p® p® 


FIGURE 13.0 JOINT INTERNAL FORCE SYSTEMS 


reliance is placed on teat observations to formu- 
late the approsiinate lineax theory presented. 

The failxire ^de of typical fastener critical 
single shear tests of high strength titanium 
fasteners in an aluminum Joint is nearly always 
tension, the failure occurring at the root area of 
the fastener or by stripping the nut. Also, post 
inspection of the failed fastener usually reveals 
no evidence to suggest excessive bending or shear 
deformation. This behavior indicates the tension 
load is high and is of sufficient magnitude to 
keep the fastener essentially straight. Because 
of this, the tension term is assumed to be the 
predominate term in the fastener stiffhess equa- 
tion. The bending and shear terms, being small, 
are neglected and the fastener assumed rigid for 
these modes of deformation. 

figures 12 and 13 define the geometry and 
internal force systems considered in the deriva- 
tion of the tilt factors. The assumption that the 
fastener is rigid for bending and shear permits a 
linearly varying bearing stress distribution to be 
used. 


The resulting equations for the tilt factors 
are as follows; In the interest of brevity, only 
the final eqiiations used to calculate the tilt 
factors are presented. By definition 


e 


Peak Hole Bearing Stress 
Uniform Hole Bearing Stress 



br 


for element ( 1 ), thicioiess (t^) 


e. 


. “ ^ * 5 S^t , (t^ . 

(x t 2) 

or element ( 2 ), thickness (t^) 

. W 3 s/2 (t, . t^)/ (^3^ , 
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7T 




(t, ♦ tj) 


iZ * -55 L, 






Where 




'•f '=•!» 


Fhstener Shank Diameter (ins) 
Fastener Head Diameter (ins) 
Thread Pitch Diameter (ins) 

Nut Outside Diameter (ins) 

D. or D , whichever is the 
^ smaller (ins) 

£^aged Thread Length (ins) 

Modulus of elasticity for Plate, 
Fastener and Nut, respectively 
(psi) 


Figure i4.0 shows the result of applying the 
tilt factor equations to a single lap shear alumi- 
num Joint with titanium TaperLok fasteners (shear 
type) and aluminum torque-off collars. The tilt 
factors 6-1 and $2 versus z /D are presented for 
t^/t^ » .5 and 1.0, ' 

Co unte rsxm k Fastener 

T?o account for the influence of a countersunk 
fastener it is assumed that the depth of tne 
countersunk nead is effective. 

The thicioiess of the countersunk plate is, 
therefore, reduced to: 


« t - D 


esk/_ 


The reduced thickness (t' ) is then used in 
place of the actual thickness (t), to enter 
Figure l4.0 for 6 and to calculate the uniform 
hole bearing stress 


In the above equations is the relative 
stiffhess of the plate (k ) to the total stiffness 
of the plate and fastener^(k -► k^) 







and k- 


Cb- Ej 




The term C is a pre-load factor to account for 
fastener pretension load caused by nut installation 
torque and b is the fastener effective moment arm 
defined by the equation 


b 


Ji 

37T 


- p^) 

(P32. p 2 ) 


The term A is the equivalent fastener shank 
area over length (t Z-) and accounts for fas- 
tener thread and nut effects 


5.0 Test Correlation 

To correlate the fatigue analysis method 
presented and test, the results from two full 
scale fatigue test articles, X99? and X99b, are 
used. The strucnire considered is the inner wing 
box lower surface assembly between the wing root 
rib and the inboard engine pylon. 



FIGURE U.O tilt factor SINGLE LAP SHEAR 

ALUMINUM JOINT WITH PROTRUDING 
HEAD titanium TAPER-LOK FASTENERS 
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The dpectrua of loads applied to the Z993 A&d 
X996 test articles were based on the block method 
of fatigue testing. The X993 an expedited 

wing fatigue test -and designed to yield early data 
on potential wing fatigue problem areas. Uie 
block spectrum of loads was, therefore, reduced to 
a ns miTn and had ^ load coxuiitions with a total 
of i6l ,000 cycles per lifetime of >0,000 cyclic 
test hours. The X99& test was ^e primary wing/ 
fuselage fatigue test and had 16 load conditions 
with a total of ^00,000 cycles per lifetime of 
>0,000 cyclic test hours. In both the X99> and 
X99S test spectra the loads were applied in passes 
of >000 cyclic test hours, and the load sequence 
was low-high- low. Also, using the K-,^thod of 
fatigue analysis, with notched 3-K data, the X99> 
and X99o test spectra were derived to give the 
same analytical damage for a * ^.0. 

Figure 1>.0 shows a planform view looking 
down on the inner wing box lower surface assembly 
and indicates the specific analysis points used in 
the correlation. The assembly shown consists of: 

3 spar caps, 9 surface panels and the attachment 
of 3 spar webs. The assembly contains 1> spanwise 
single lap shear joints. The spar caps and surface 
panels are machined from 7C75-T63V1 extrusion, and 
the spar webs are 7^75-T65 i plate. The spar cap-to- 
panel and panel -to -panel joints are formed using 
countersunk head titanium TaperLok fasteners 
with aluminum collars. The spar web-to-apar cap 
joints are formed using protruding head titanium 
TaperLok fasteners with alimanun collars. IXiring 
assembly the faying surfaces of all joints are 
coated with sealant and ail fasteners are wet 
installed. 

Table 1 .C lists, for the 5 analysis points 
shown on Figure 15.0, the joint geometry and the 
empirically derived factors K^, K. and a ^ . The 
empirical factors are based oh xhi teat demon- 
strated fatigue life at each point and are calcu- 
lated for X993 and X99a using the following 
methods: 

a), K,. - Conventional K_-method of fatigue 
analysis using^7-73-T6 notched 3-K 
data and considering uniaxial stress 
only. 


IWBR STA IWBR STA IW6R STA 

M6.B0) 246.104 366.191 



FIGURE 15,0 PLANFORM VIEW OF INNER WING BOX 
LOWER SURFACE PANEL ASSEMBLY^ CORRELATION' 
ANALYSIS POINTS 


5), K. - Conventional K,^-method of fatigue 
analysis using* the 7073-T6 joint 
3-:; data piresented in Section 2.j 
and considering uniaxial stress 
only . 

C),a3 • SSF-roethod of fatigue analysis pre- 
sented in Section 2.C, considering 
uniaxial stress and shear flow or 
load transfer. 

The empirical factors Y^. and K. are the same 
as the ’quality level index - X* (^^*^with the sub- 
scripts K and J added to denote the type of S-IC 
data used in their derivation. The empirical 
factor fl/3 is the combined a and 0 factors dis- 
cussed in Section 2.0'. On the full scale fatigue 
test articles it was impossible to separate the 
effects of hole preparation and hole filling; 
therefore, the a and 0 factors are combined and 
referred to as a /3 . 

Comparing the empirical factors for X993 and 
199 B shows that the req\iired correlati:n factors 
(X99B/X993) are a maxlrmun for method A and a 
mi n’iffmm for method C, the required correlation 
factors for method 3 falling between these extreme 
values. For example, at point (^) the correlation 
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METHOD A 
METHOD B 


□ 


method c 



ON X993 EMPIRICAL FACTOR, TO 

demonstrated life 


factors are T.2 j, 1.1^ and 1.07 for methods A, B 
and C. respectively. If the average values are 
used this trend is maintained » the correlation 
factors being i-i5. 1.09 and 1.03- 

Tp further compare the empirical factors the 
scatter about the overall average values shows 
that the scatter xising method A is •*•2358 to -1358 
of the overall average values of 5*02, the scatter 
using method 3 is to -1058 of the overall 

average value of 4.21 and the scatter using method 
C is •*•9)8 to -S5E of the overall average value of 
-. 9 ^ 5 . 

A comparison of the empirical factors does not 
necessarily indicate the true accuracy of a method 
of fatigue analysis « as a small change in the 
empirical factor can produce a large change in 
predicted life. Figure 16.O illustrates this 
point and shows the ratio of the X99S predicted 
life, based on the X993 empirical factors, to 
demonstrated life. I3ie trend is the same as the 
trend observed for the empirical factors, the life 
predictions being in the range 1.6 to 2.3» 1*2 to 
1.3 and 0.9 to 1.3 times the demonstrated life for 
methods A, B and C, respectively. Ibe best life 
prediction for method B of 1.2 and the worst life 
prediction for method C of 1.3 occur at point (2); 
at all other points method C gives the best life 
prediction. 


6.0 Conclusions 


The correlation of analysis and test results 
presented in Section shows that the beat cor- 
relation is achieved using method C. As method 
C is the SSF-method of fatigue analysis and in- 
cludes fastener load transfer effects, the factor 
used to define the fatigue quality of a Joint is 
not only a function of the Joint geomatry, but 
also of the Joint loading. This conclusion then 
points to the need axid iagsoxtanoe of such an 
analysis which considers, say on a %iing box, not 
only the stress due to wing bending, but also the 
effects on fatigue life of wing shear and torsion. 
This method then should be considered a valuable 
tool in developing fatigue test spectra and when 
comparing fatigue design spectra on structure vrith 
different usages. On the C-3A wing, the SSF-method 
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of fatigue analysis, for single lap shear Joints, 
is being used to determine the effects on the wing 
structure of different methods of lift distribu- 
tion control systems and/or fuel utilization 
sequencing. 
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SOFIA 



CRITERIA AND ASSUMPTIONS 

• ALLOWABLE STRESS LEVELS FATIGUE CRITICAL 
ALUMINUM PARTS 

- FLIGHT & PRESSURIZATION LOADS - 45KSI ULT 

- FLIGHT & PRESSURIZATION LOADS UNINSPECTABLE AREAS - 
30KSI ULT 

- PRESSURIZATION LOADS ONLY -21 KSI ULT 

• ABOVE GROSS AREA STRESSES BASED ON 

- SECTION EFFICIENCY 75% 

- FATIGUE QUALITY INDEX 4.5 

• FOR OUTSTANDING FLANGES WITH NO HOLES 

- PRESSURIZATION LOAD ONLY - 45KSI ULT 















SOFIA 

CRITERIA AND ASSUMPTIONS 

* 

• ALLOWABLE STRESS LEVELS 

- NETSECTION ALLOW YIELD AT LIMIT LOAD 

- NET SECTION ALLOW ULTIMATE AT ULTIMATE LOAD 

• STIFFNESS 

- MAINTAIN SAME STIFFNESS AS BASIC 747 AS FAR AS FEASIBLE 

• TELESCOPE MASS PROPERTIES 

- USE LMSC VALUES FROM LMSC/P01 6023 


BASELINE FEA MODEL 

• DISCRETELY MODELED FS1480-FS21 80 

- 5000 NODES 

- 31000 DOF 

• SUBSTRUCTURE FWD OF FS1480 

- REACTIONS AT STUB WING TIPS AND FORWARD END OF 
FUSELAGE (FS900) 
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SOFIA 


I El G e 1 n c| 7 7 C e n t. e r -Se c t. i o n 5 li 1'j s t r' u c t u. r e 
IFS 90lZ) — F5 1480 Including Stub Wings 
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SOFIA FEA MODEL 
• BASELINE 747 MODEL PLUS 

- FORWARD BULKHEAD ALL ELEMENTS DISCRETELY MODELED 

- AFT BULKHEAD BEAM ELEMENTS AND SHEAR WEBS 

- TRANSITION DISCRETELY MODELED DOWN TO FASTENER 
LEVEL AT ONE LOCATION - OTHER LOCATIONS MODELED AS 
BEAMS 

- SURROUND STRUCTUREiBEAMS, RODS & SHELL ELEMENTS 

- FLOORING BEAMS, RODS AND SHEAR ELEMENTS 



LOADS 

• DEVELOPED FROM ENVELOPE LOADS IN BOEING 
REPORTS 

• LOAD CASES 

- VERTICAL BENDING ENVELOPE 

- LATERAL BENDING ENVELOPE 

- PRESSURE ONLY 

- LATERAL LOAD - LEVEL FLIGHT 

- 2P CONDITION ‘ 

- BENDING CASES PLUS ULTIMATE PRESSURE 





Shear Flow (Ib/in) 



LOADS/MODEL VALIDATION 

Aft fuselage without Doors FS 1740 Aft Fuselage without Doors FS 1740 




747-100 \7eight Conditions 


Cond . 

GW 


c.g. 


Izz 


V'.ips 


MAC 

in“lb-sec2 

in-lb-cec2 

000 

327 


34 .9% 

.250E+09 

.374E+09 

001 

360 


13.0% 

.305E+09 

.422E+09 

002 

585 


13.0 

.480E+09 

.611E+09 


747-100 LIMIT 

Acceleration 

Designt Conditions 

(UV, 

M'ddot apply to Lateral Gust, Hz, Oddot to others, 

Gust Conditions at Vc; 

375 

KEAS at 

20kFT, or HO. 92 

at 40 kFT] 

Cond . 

Type 

GV/ 

c.g. 

Alt. Hy/Hz 

u/'l' 



kips 

IMAC 

KFT LIMIT 

rad/scc2 

GLOl.l 
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0.13 
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0.100 

GNOl , 1 
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0.13 

40 -0.16 

0.201 
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20 3.02 
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0.13 
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2WLL 1 
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0.13 

S.L. ' 2.15 

-0.190 


2 WLL Two-wheel Level Landing at Maximum Vertical Reaction 


MAC 


MAC 

327.8 in 



LEMAC 

.1258 in 



S 

5500 ft-sq 
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. 2348 in 



LFUS 

2702 in 
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2495 in 



Ltai 1 

1254.05 in 
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1470.7 ft-sq 
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6,96 
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; 3.6 
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. 20 
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Ve 

375 
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CLaC to 
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1/deg 

Cmat to 
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' -1.0258 

1/rad 
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0.1964 

0.1964 

1/rad 



FORWARD BULKHEAD 

* 

# DISCRETELY MODELED 

- BULKHEAD WEB 

- BEAM WEBS 

- BEAM CAPS 

- BULKHEAD CAP 



• • • 

SOFIA f\JASA 


FORWARD BULKHEAD 

• DEFLECTION 0.5” LIMIT AT TELESCOPE MOUNTS 

• BEAMS 20 INCH DEEP 

• BEAM WEBS NON-BUCKLING AT LIMIT LOAD 

• EXTRUDED PRESSURE PANEL 

• BEAM CAPS 5.6 IN2 

• BEAM WEBS 0.063 AND 0.080 


••es 



SOFIA 


SCFIfH - f\ixi bulkhead beam cap axial stres; 
forujard caps (FS I7E0) 

I P C9 ps i ) Cond i t ion 





90FIR - fiajd bulkhead beam uieb shear floes 
( Ib/'in) 

GL0S. 1 Lateral Guet 0 40k ft ^ brake off 
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SOFIA 



AFT BULKHEAD 

• 10 INCH DEEP BEAMS 

• EXTRUDED PRESSURE PANEL 

• 0.063 BEAM WEB 

• 2.0 IN^ BEAM CAP 

• STIFFNESS COMPATIBLE WITH FORWARD 
BULKHEAD 
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BULKHEAD TO SHELL TRANSiTION 

• AREAS AND INTERCOSTAL DEPTH VARY FOR 
SMOOTH TRANSITION 

• LOW STRESS LEVELS 

• LOW LOAD TRANSFER RATE 

• LOW STRESS SEVERITY FACTOR 

• LONG LIFE 
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Elem Area Load Stress Fastener 
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CUT-OUT SURROUND STRUCTURE 

• CONVENTIONAL CONSTRUCTION 

• PRIMARY (EXTERNAL) AND SECONDARY 
(INTERNAL) LONGERONS 

• DOUBLED UP SKIN 

• TRIPLER FOR ADDED REINFORCEMENT AT 
CUT-OUT EDGE 

• LONGERONS SPLICED WITH TENSION FITTINGS 
AT BULKHEAD 

• STRINGERS SPLICED WITH SHEAR JOINT AT 
BULKHEADS 
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FLOORING 


• MAIN DECK FLOOR REMOVED IN CAVITY 

• MAIN DECK & CARGO DECK MAINTENANCE 
WALKWAY 

• CUT-OUT FORWARD OF BULKHEAD 

• FLOORING ACTS AS BEAM TO HOLD SHELL 
SHAPE 

• FLOORING DIVORCED FROM BULKHEAD TO 
PERMIT BULKHEAD DEFLECTION 


SC-r L-* “ floor forward of f ud bulkheao 
otfleotionc 5Qx / axial forces 
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FUTURE TASKS 

• CONTINUE DEVELOPMENT OF DESIGN IDEAS 
DEVISED DURING PHASE I & II 

• REQUIRES MORE DETAILED DATA ON 

- TELESCOPE SUPPORT 

- PRESSURE SEAL 

- TELESCOPE BRAKE 

- ORIFICE DOOR & MECHANISM 
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GOST IMPACT OF VEHICLE SELECTION 

• AIRCRAFT COST 

- 747-SP $10 - $12 MILLION 

- 747-100 $2.5 - $4 MILLION 

• LOADS AND MODEL DEVELOPMENT 

- 747-SP NO DATA $5 - $10 MILLION 

- 747-SP FUSELAGE STRESS REPORT $2 MILLION 

- 747-SP FUSELAGE STRESS & LOADS REPORT $1 MILLION 

- 747-100 $300,000 



